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I .  Introduction 

During  the  past  year,  our  research  in  "white-light  otpical  information 
processing  and  holography"  has  demonstrated  several  significant  results 
that  several  of  optical  information  operations  can  be  carried  out  with  non¬ 
coherent  source.  As  in  the  past  year,  we  have  been  quite  consistent  in 
reporting  our  AFOSR  sponsored  research  in  various  journals  and  conference 
publications.  Sample  copies  of  these  papers  are  included  in  this  report 
as  in  the  following  chapters  to  provide  concise  documentation  of  our  work. 
In  the  following  sectins ,  we  will  provide  a  general  overview  of  our 
research  progress  made  in  the  past  twelve  months.  Details  on  some  of 
those  research  progress  are  immediately  provided  in  subsequent  sections. 

A  list  of  publications  resulting  from  AFOSR  support  is  included  at  the 
end  of  this  report. 

II.  Summary  and  Overview 

Since  the  intended  research  program  is  for  the  white-light  optical 
information  processing  and  holography,  attention  gathering  the  optical 
components  for  use  in  such  a  processing  system  is  necessary.  The  basic 
system  involves  white-light  (Xenon  and  Zirconium  arc  lamps),  incoherent 
(Mercury  arc  lamp)  light  sources,  achromatic  lenses,  diffraction  gratings, 
(amplitude  or  phase  type) ,  complex  spatial  filters  and  various  supporting 
electro-optical  elements.  The  basic  systems  of  this  white-light  or 
incoherent  light  optical  information  processing  system  was  synthesized 
during  the  earlier  part  of  the  first  six  month  period,  as  shown  in  Fig.  1. 
Since  the  use  of  incoherent  source  for  the  processing,  the  physical 
processing  environment  is  not  ns  critical  ns  coherent  optical  processor: 

For  example,  weighty  optical,  bench  Is  not  required  and  it  is  not  required 
to  have  completely  dust  free  environment.  This  white-light  optical 
processor  greatly  simplifies  the  stability  and  ruggedness  of  the  pro- 
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cessing  system.  Another  interesting  feature  of  this  system  is  that  the 
liquid  gates  are  generally  not  required  for  the  processing  operation. 
Futhermore,  because  the  use  of  broad  spectral  band  white-light  source, 
the  system  is  particularly  suitable  for  color  image  processing. 

2.1  Smeared  Image  Deblurring  (See  Section  III) 

Our  research  in  the  past  year  period  has  included  the  smearud- 
photographic-image  deblurring  [1],  sample  of  the  result  is  shown  in 
Fig.  2a.  In  comparison  with  the  result  obtained  by  coherent  source, 
the  one  obtained  with  the  white-light  source  offers  a  fewer  artifact 
noise.  I.  wever,  the  deblurring  spatial  filter  that  we  used  were  a 

o 

narrow  spectral  band  centered  at  5154A  green  light.  To  compensate 
for  the  scaling  of  the  signal  spectrum  due  to  the  wavelength  of  the 
white-light  source,  a  fan-shaped  deblurring  filter  should  be  used. 

The  generation  of  a  fan  shape  deblurring  filter  by  means  of  computer 
technique  and  optical  coating  technique  will  be  reported  in  the 
future  annual  reports. 

We  have  also  extended  this  Image  deblurring  technique  to  color 
photographic  image  deblurring  [2],  as  shown  in  Fig.  3.  We  have  pointed 
out  earlier,  the  white-light  processing  technique  is  suitable  for  color 
image  processing,  the  color  image  deblurring  was  achieved  with  narrow 
red  and  green  color  deblurring  filters.  The  center  wavelengths  of  this 

O  O 

set  of  deblurring  filters  were  about  6328A  and  5461A  respectively.  The 
filter  bandwidth  contained  five  main  lobes  and  thus  spectral  widths 

O 

are  about  100A.  From  the  result  of  Fig.  3,  we  see  that  the  deblurring 
effect  of  the  green  letters  "P"  and  "C"  seems  to  be  more  effective  than 
other  color  letters.  This  is  primarily  due  to  linear  amplitude  transmittance 
of  these  smeared  letters.  In  other  words,  the  smeared  letters  of  red  and 
yellow  colors  were  more  saturated. 


Figure  2.  Linear-Photographic  Image  Deblurring 


Figure  3.  Color  Image  Deblurring 


We  note  that,  the  results  of  Image  deblurring  were  utilized  narrow 
spatial  filters  in  the  spatial  frequency  plane.  Extension  to  the  whole 
spectral  band  of  the  white-light  source  on  the  aspect  of  the  synthesis 
of  a  fan-shape  spatial  filter  concept  for  image  deblurring  is  the 
subject  of  current  research  that  will  be  reported  upon  in  the  future 
annual  reports. 

2.2  Image  Subtraction  (See  Section  IV) 

We  have  in  the  past  year  investigated  the  possibility  of  using 
incoherent  and  white-light  sources  for  image  subtraction.  Since 
optical  image  subtraction  is  a  one-dimensional  processing  operation, 
instead  of  using  a  point  source  of  light,  a  line  source  can  be 
utilized.  Moreover,  the  image  subtraction  operates  upon  the  one-to- 
one  correspondent  Image  points,  a  strictly  broad  coherence  requirement 
is  not  needed.  It  is  possible  however  to  encode  an  extended  incoherence 
source  to  obtain  a  point-pair  spatial  coherence  function  for  the  image 
subtraction  operation  [1-5]  (sec  section  V).  In  evaluating  the  spatial 
coherence  requirement,  we  had  applied  the  partially  coherent  theory. 

In  experimental  demonstration,  we  provide  a  set  of  continuous  tone 
image  transparencies  as  input  objects  as  shown  in  Figs.  Aa  and  Ab .  By 
comparing  these  two  figures,  we  see  that  a  liquid  gate  was  withdrawn 
from  the  optical  bench  in  Fig.  Ab .  Figure  Ac  shows  the  subtracted 
image  obtained  with  the  incoherent  source  encoding  technique,  while 
Fig.  Ad  was  obtained  with  coherent  processing  technique.  From  Fig.  Ac, 
a  profile  of  subtracted  liquid  gate  can  readily  be  seen,  which  from 
the  result  obtained  with  coherent  source,  the  subtracted  image  is 
severely  damaged  by  the  coherent  artifact  noise.  Thus  we  see  that  the 
incoherent  technique  is  indeed  offering  a  better  image  quality. 


Obtained  With  Incoherent  Technique. 
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We  have  also  extended  this  image  subtraction  technique  for  color 
images  [6,7],  as  shown  in  Fig.  5.  Fig.  5a  and  5b  shows  two  color 
images  of  a  parking  lot  on  input  color  transparencies.  We  note 
that  a  red  color  passenger  car  in  Fig.  5a  is  missing  in  Fig.  5b. 

Figure  5c  is  the  color  subtracted  image  obtained  with  the  incoherent 
source  encoding  technique.  In  this  figure,  a  profile  of  red  passenger 
car  can  clearly  be  seen  at  the  output  image  plane.  It  is  also 
interesting  to  see  that  the  parking  line  (in  yellow  color)  on  the 
right  side  of  the  red  car  can  readily  be  seen  with  the  subtracted 
image . 

Again,  we  note  that,  the  results  of  image  subtraction  with  encoded 
extended  incoherent  source  were  accomplished  by  narrow  spectral  band 
of  light  sources.  Extensive  additional  research  on  the  aspect  of 
source  encoding  for  extended  white-light  source  is  the  subject  of 
current  research  that  will  be  reported  upon  in  future  annual  reports. 


2.3  Visualization  of  Phase  Object  (See  Reference’  8) 

We  have  extended  the  color  image  subtraction  technique  for  visualiza¬ 
tion  [8]  as  shown  in  Fig.  6.  By  this  variation  of  color  coded  fringe 
pattern,  it  is  possible  to  determine  the  phase  variation  in  more  detail. 


The  detail  analysis  of  this  work  will  be  reported  in  our  next  annual 


report . 


2.4  Coherence  Requirement  (See  Reference  9) 

We  have  also  in  this  period  evaluated  the  coherence  requirement 
for  white-light  optical  processing  [9].  The  spatial  and  temporal 
coherence  requirements  for  smeared  Image  deblurrlng  and  image  sub¬ 
traction  problems  were  evaluated.  For  image  deblurring,  we  show 
that  spatial  coherence  requirement  Is  dependent  upon  the  source  slit 
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Figure 


jn"J 
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size,  the  smeared  length,  and  the  spatial  frequency  of  the  grating. 

For  image  subtraction,  the  spatial  coherence  requirement  is  dependent 
upon  the  ratio  of  slit  size  and  the  spacing  of  the  slit  of  a  source 
encoding  mask,  spectral  bandwidth,  and  the  separation  of  the  input 
objects.  From  the  obtained  coherence  requirements,  we  were  able  to 
develop  a  general  source  encoding  concept  for  partially  coherent 
optical  processing  [3]. 

The  detail  analysis  of  the  coherence  requirement  will  be  reported 
in  our  next  annual  report. 

2.5  Remarks 

We  have  demonstrated  experimentally  white-light  optical  informa¬ 
tion  processing  technique  utilizing  a  diffraction  grating  method.  As 
compared  with  the  conventional  coherent  optical  processing  technique, 
the  white-light  processing  technique  offers  several  basic  advantages: 

1.  The  cost  of  the  processing  system  is  significantly  lowered 
with  the  elimination  of  the  laser  sources. 

2.  The  alignment  procedure  of  the  white-light  processing 
system  is  generally  simpler  than  the  coherent  technique. 

3.  The  output  results  are  free  from  artifact  noise  that 
generally  plagues  the  coherent  processing  systems. 

In  short,  the  white-light  optical  information  processing  system 
is  simple,  versatile  and  economical  to  operate.  We  would  emphasize 
that,  there  are  many  information  processing  operations  which  can  be 
carried  out  rather  easily  with  a  white-light  source.  In  other  words, 
if  a  hi gh  spatial  coherent  requirement  (a  requirement  lias  to  be 
evaluated)  for  an  information  processing  operation  is  not  needed,  then 
a  white-light  processing  technique  will  be  more  advantageous.  However, 
if  the  demand  of  the  spatial  coherence  for  an  information  processing 
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operation  is  high,  then  coherent  sources  may  not  be  avoided.  We  note 
that,  the  analysis  of  the  cost  and  merit  as  a  function  of  coherent 
requirement  for  the  proposed  white-light  processing  technique  will  be 
evaluated  in  our  future  research. 

2.6  Future  Research 

Our  future  research  is  expected  to  follow  the  general  direction 
addressed  in  the  past  year  with  a  major  attention  given  in  the  following: 

1.  We  will  synthesize  a  matched  filter  that  is  suitable  for  the 
white-light  optical  processing  technique.  Experimental 
demonstrations  for  the  incoherent  and  the  white-light  processing 
technique  will  be  carried  out. 

2.  We  will  carry  out  a  computer  generated  spatial  filter  program 
that  is  suitable  for  the  broad  spectral  band  optical  processing 
program.  In  other  words,  a  multiwavelength  and  multiband 
spatial  filter  may  not  be  generated  by  coherence  sources. 
However,  to  compensate  the  wavelength  variation  it  is,  in 
principle,  to  be  generated  by  computer  techniques. 

3.  We  will  carry  out  the  experimental  confirmations  for  those 
computer  generated  spatial  filters  for  complex  signal 
detection,  image  deblurring,  image  subtraction,  color 
encoding  problems,  etc.  for  a  broad  band  of  white-light 
source . 

4.  We  will  investigate  a  source  encoding  technique  for  the  image 
subtraction  for  broadband  white-light  source.  Experimental 
confirmations  will  concurrently  carry  out  for  monochrome  and 
color  image  subtraction. 

5.  Since  the  white-light  source  contains  all  the  visible  wave¬ 
length,  we  will  develop  a  research  program  primarily  devoted 


13 


to  color  Image  processing.  We  note  that,  initially  all 
the  visible  images  are  color. 

6.  We  will  investigate  the  possibility  of  implementing  natural 
sun  light  in  the  white-light  optical  processor.  Experiments 
with  the  sun  light  processing  technique  will  concurrent 
carry  out.  We  expect  interesting  results  will  be  obtained 
with  the  sun  light  source. 

7.  The  analysis  of  the  cost,  merits,  and  limitations  of  the 
proposed  white-light  processing  technique  will  be  evaluated. 

With  reference  to  the  analysis,  more  sophisticated  processing 
operation,  as  applied  to  wide-band  signal  and  digital  signal 
processings  will  be  carried  out. 

8.  We  will  investigate  the  possibility  of  developing  a  real-time 
white-light  optical  processor.  This  real-time  processor  should 
have  the  capability  of  performing  the  major  optical  processing 
operations  that  a  coherent  processor  can  provide.  We  will  also 
investigate  the  possible  application  of  the  white-light  and  sun 
light  processing  technique  to  some  problems  in  optical  computing 
(e.g.,  digital  or  analog). 

9.  We  would  also  investigate  various  possible  applications  of  the 
white-light  optical  signal  processing  for  tracking  and  identifi¬ 
cation  of  aircrafts  and  missiles;  spread  spectrum  communications; 
counterfeit  deterrence;  acoutical-optical  sonic  spectral  analysis, 
and  recognition;  aerial  and  satellite  picture  processing, 
enhancement  and  identification;  and  application  to  various  bio¬ 
logical  and  medical  aspects. 

In  short,  our  goal  is  to  develop  a  full  research  program  on  the 
white-light  signal  processing  technique  which  has  the  capability  of 
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carrying  out  major  processing  operations  that  a  coherent  optical 
processor  can  offer.  This  white-light  optical  processing  system 
is  potentially  important  in  many  areas  of  monochrome  and  color 
signal  processings.  We  stress  that  this  white-light  processing 
technique  will  provide  a  first  step  toward  the  research  and 
development  of  white-light  and  sun-light  optical  computers.  We 
note  that,  if  the  demand  of  the  spatial  coherent  requirement  for 
an  information  processing  operation  is  not  stringent  (which  will 
be  determined  in  our  future  research  program) ,  then  the  white-light 
signal  processing  technique  is  more  convenient  and  reliable  than 
the  conventional  coherent  optical  processing  techniques.  However, 
if  the  spatial  coherent  requirement  becomes  very  stringent,  then 
the  use  of  a  coherent  source  may  not  be  avoided. 
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Abstract.  Recent  advances  in  electro- optics  have  brought  into  use  communica¬ 
tion  and  information  theory  to  analyze  performance  in  coherent  and  incoherent 
optical  information  processing  systems  An  optical  information  processsing 
system  can  be  analyzed  with  many  of  the  same  concepts  of  linear  system  theory 
(eg,  spatial  impulse  response,  spatial  frequency  and  spatial  domain  synthesis, 
etc  ),  and  the  photographic  images  to  be  processed  can  be  regarded  in  the  same 
manner  as  time  signals  (e  g  spatial  frequency  content,  spatial  amplitude  and 
phase  modulation,  space-bandwidth  product,  etc.)  Both  coherent  and  incoher¬ 
ent  optical  processing  systems  can  be  treated  as  linear  systems,  and  the 
processing  operation  can  generally  be  carried  out  by  communication  theory 
concepts 

Although  coherent  optical  information  processing  operations  have  been 
used  for  performing  complex  amplitude  operations,  complex  processing  can 
also  be  performed  with  incoherent  or  white-light  illumination  The  importance 
of  optical  information  processing  operations,  either  coherent  or  incoherent,  is 
due  to  the  basic  Fourier  transform  properties  of  lenses  In  this  paper,  we  will 
discuss  mostly  the  incoherent  systems  because  they  are  of  more  recent  interest 
and  possess  certain  advantages,  we  feel,  over  the  traditional  coherent  optical 
processors  Experimental  illustrations  of  the  results  are  provided 

In  view  of  the  broad  area  in  optical  processing  of  photographic  images,  we 
will  confine  ourselves  to  a  few  applications  that  we  consider  of  general  interest. 
We  apologize  for  the  omission  of  other  techniques  and  applications,  and  for 
neglecting  the  inclusion  of  their  references 
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1.  INTRODUCTION 

Communication  and  information  theory  was  originated  by  a  group 
of  mathematically  oriented  electrical  engineers  whose  interest  was 
centered  on  electrical  communication  Nevertheless,  from  the  very 
beginning  ol  this  discos ei  >  ol  communication  and  mloimalion  the¬ 
ory.  interest  m  its  application  to  optical  systems  has  been  sigotous 
As  a  result  ot  i event  advances  m  optical  signal  pioccssing  and  optical 
communications,  the  relationship  between  optics  and  communica¬ 
tion  theory  has  grown  very  rapidly. 

Mention  must  he  ntadeol  a  lew  impoilant  early  coiiltihiilions  to 
this  field  It  was  in  the  early  1950s  that  the  communication  and 
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information  theory  aspects  of  optical  processing  techniques  first 
became  ev  idem  I  lie  most  important  impact  must  be  due  to  (labor's 
work  on  light  and  infoi illation1  m  1951. 1  has.  (irey.  and  Robinson’s 
woik  on  I  mu  ici  tiealment  ofoptic.il  processes-  in  1952.  Elias's 
pa  pc  i  on  optics  and  communication  theory  '  in  1955;  and  I  ornldo  di 
I  rancia's  woik  on  tcsolwng  power  and  information4  in  1955  How¬ 
ever.  the  veiv  lust  application  of  communication  theory  to  modern 
optical  mini matioii  pioccssing  was  piohahly  O'Neill's  woik  on  spa¬ 
tial  liltcimg  in  optics'  in  195b.  Recause  ol  the  broad  mtcicsl  in  this 
field  at  that  lime,  a  special  symposium,  ('ommunication  and  Infor¬ 
mation  I  lieoiy  Aspect  ol  Modern  Optics,"  took  place  in  I960  Since 
then,  the  application  of  communication  and  information  theory  to 
optical  ignal  processing  has  commanded  great  inlctcst  I  he  applica¬ 
tions  ot  cotic.d  spatial  filtering  were  particularly  evident  in  the  field 
of  radar  sig  iuI  processing,  and  it  was  in  this  field  that  Cutrona.  Leith. 
Palermo,  and  Porcello  published  a  classic  article  on  optical  data 
processing  and  liltcimg  sy  terns’  in  I960  fins  iitide  stimulated  a 
bioail  mtcicsl  hi  optical  processing  ol  pliotogiapliic  images  With 
the  invention  ol  a  stiong  coherent  source.  i.e„  laser,  in  the  early 
I9h(h,  I  citli  and  llpalnieks'  woik  on  icconsliuctcd  w-aveliout  und 
communication  theory"  allowed  lor  the  first  time  the  formation  of 
high  quality  holographic  images.  I  Ismg  the  spatial  frequency  carrier 
concept  ol  I  villi  and  Upalnicks  lot  hologiaphy.  the  I9h4  paper 
‘’Signal  Selection  by  Complex  Spatial  Filtering"  by  Vandcrl.ugt* 
introduced  the  subject  of  optical  character  recognition  via  the  optical 
matched  filter  correlator  Since  then  optical  information  processing 
has  been  applied  to  a  large  variety  of  problems  Thcrefoie.  it  is 
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evident  that  communication  and  mtnrmalion  thcorv  has  stimulated 
a  bio. id  range  ul  application  In  moduli  optical  mini mai ion 
processing 

f'oheiem  optical  mini  illation  processing  opeialmns  have  tradi¬ 
tionally  been  legatdcd  av  the  most  useful  lor  palm ming  complex 
amplitude  upciulions  llmvcsci,  ill  icccnt  vc.ux  complex  inloima- 
lion  processing  performed  with  mcohcient  or  white  light  illiimina- 
tion  lux  iceeived  meieaxed  attention  Ncveithclexx.  the  inipoitanee 
ol  optical  mini mation  piocessing  opeiationx.  eitliei  coherent  or 
incohetent,  is  pi  imarily  based  on  the  basic  I  miner  transform  proper- 
tv  ol  lenses  l  he  applications  ot  these  two  optical  processing  systems 
to  phoioginphic  images  will  be  described  I  x  pci  internal  illusti.it  ions 
ot  the  results  will  he  provided 

In  view  ol  the  hioad  area  ill  optical  pioccxxmg  ol  phoioginphic 
images  and  the  coniines  ol  space,  we  will  discuss  a  lew  applications 
that  w e  consider  ol  general  interest.  I  he  field  has  now  matured  to  the 
point  where  there  are  mans  good  textbooks  available,  and  we  refer 
the  general  reader  to  seveial  of  these.1"  " 


this  section  we  will  describe  a  white-light  processing  technique  by 
which  phoioginphic  images  can  be  piocexxcd  with  complex  ampli¬ 
tude  tillci  s  Wc  iioic  that  this  w  lute- light  piocexxor  may  a  He  v  talc  the 
m.iioi  disadvantages  imposed  on  the  coherent  processing  system. 

Now  we  desenbe  a  photographic  linage  processing  technique  that 
can  lie  can  led  out  hv  aw  lute- light  source.  4  as  show  n  in  I  ig  I.  We 
note  that  the  while-light  processing  system  is  similar  to  that  ol  a 
cohcicnl  piocexsmg  system  except  that  a  white-light  souice  and  a 
high  ditliaclion  elticiency  grating  are  inserted  m  the  input  plane  I*, 
It  we  place  a  photographic  transparency  s(x.y)  in  contact  with  the 
diffraction  gtatmg.  then  the  complex  light  field  for  every  wavelength 
A  behind  the  transloim  lens  1  ^  is 


I  (p.q. A) 


+  costp^xlj 


exp[  i(p.x  -t  qy)]dxdy  . 


(I) 


2.  INCOHKRKNT  OPT1CA1.  PROC  FSSINf.  SVSTKMS 

I  he  use  ol  coherent  light  enables  optical  xy  stems  to  c.ii i y  out  many 
sophisticated  information  processing  operations. 14  However,  coher¬ 
ent  optical  processing  systems  are  plagued  with  coherent  artifact 
noise,  which  frequently  limits  the  processing  capability  of  the  sys¬ 
tems  Although  many  optical  information  processing  operations  can 
be  implemented  by  systems  that  use  incoherent  light.1'  "■  there  are 
otliei  serious  problems.  On  one  hand,  the  incohetent  processing 
system  is  capable  ol  reducing  the  mev  it  able  attilaci  noise,  but  on  the 
other  hand  it  generally  introduces  a  dc-biax  buildup  pronlcm.  which 
again  can  lesiilt  in  pom  noise  pcilniniance  llieic  uic.  howevei. 
techniques  which  have  been  developed  loi  cohcicnl  opeialion  with 
light  ol  reduced  coherence1'  oi  noise  averaging.1' hut  at  the 
expense  ol  increased  system  complexity. 

Attempts  at  i  educing  I  lie  tc  input. d  coherence  requuetueutxoiuhe 
light  source  in  the  optical  inlormatioii  piocessing  tall  into  two  general 
categories:  (II  the  use  ol  incoherent  instead  ol  coherent  optical 
piocessing  pursued  by  I  invent  ha  I  and  Chav  cl-1  and  I  ohmann,--' 
among  others:  and  (2)  the  reduction  ol  coherence  while  still  operat¬ 
ing  with  Imcur-m-nmplitiidc  systems  ptusued  by  I  eitli  and  Roth  ■ 1 
I  lie  latter  is  the  one  that  we  believe  to  be  a  very  promising  technique, 
and  it  is  the  concept  that  \  u  is  pursuing 

Since  its  invention  the  laser  has  become  a  useltil  tool  for  many 
applications  in  coherent  optical  information  processing.  I  bis  trend 
ol  advances  was  mainly  due  to  the  complex  amplitude  processing 
capability  .  In  addition  to  the  noise  problem  mentioned  above,  the 
coherent  sources  are  usually  more  expensive  and  the  system  stability 
is  usually  critical 

Recently,  m  looking  at  optical  inhumation  piocessing  techniques 
from  a  different  viewpoint,  a  question  arose.  Is  it  necessary  that  all 
the  mini  mat  ion  pi  occxxmp  operations  k  quite  a  coho  cut  souice'’  We 
louml  that  thcie  ate  seveial  optical  mloi mation  piocessing  opcia- 
tionx  that  can  he  carried  out  with  i  educed  coherence.  II  thecoheience 
icquircmcnt  loi  certain  mini  mation  processing  npciwlious  is  not  too 
high,  it  may  be  carried  out  with  incoherent  oi  white-light  sources  In 


fig  t  A  whit#  light  optical  processor.  I.  white-light  source.  T(a).  phase 
grating:  L.  achromatic  lens:  H(p.q).  complex  spatial  filters 


where  the  integral  is  over  the  spatial  domain  of  the  input  plane  I’., 
(p.q )  denotes  the  angulai  spatial  frequency  coordinate  system,  and  C 
is  a  complex  constant. 

l  or  simplicity  of  analy  sis,  we  drop  the  proportionality  constant, 
and  I  q  (I)  becomes 

l'(p.q:AI  Slp.ql  i  S(p  -p„.q)  t  S(p  +  p(,.q)  .  (2) 

where  x(p.q)  is  the  loutier  spectrum  ol  x(x.y  I. 


I’ 


2  n 

A I 


o 


and 


v| 


2  n 
A I 


/<  ■ 


(«./l)  is  the  linear  spatial  coordinate  system  ol  (p.q)  and  I  is  the  focal 
length  ol  the  achiomatic  Hailstorm  lens  I  I  In  terms  ol  the  spatial 
coordinates  ol  o  and  /l  I  q  (2)  can  he  written. 

Af 

!  (<»./<: A>  (  ,S(ii./j)  +  C,S(«-  — 


t  C,S(||  t  p(|,/|)  . 


(3) 


liom  the  above  equation,  we  see  that  two  first -order  signal  spectral 
bauds  (i  c  .  second  and  l  hit  d  lei  ills)  ate  dispel  sed  into  lain  bow  colors 
along  the  n  axis,  and  each  spectrum  is  cetilcied  at  «  I  (Al  2n,)pll. 

In  the  analysis  we  assume  that  a  sequence  ol  complex  spatial 
lilteis  loi  v ai unis  A|(  is  available,  i  e  .  II (p„.q„f.  w licie 


I’m 


J  .1 

i 

A„l 


a  nil 


A„1 


II  we  place  these  complex  spatial  tillers  in  the  spatial  lictpiene>  plane 
uilli  each  irtiieieil  .it  o  ( A(| I  )\\r  then  the  coinplcs  lipht  lielil 
hclmul  the  spatial  licipiencv  plane  is 


I  (p.q: A)  S(p  p,,.q)  V  H(p„-  P„.qnl  (4) 

If  I 


Ol*Uf:AI  I  Nt, INI  f  III  fs|C  •  /  f'.'iitmntiHr /Oft  r  it  in  I  1  /  Vul  70  No  ?,  /  667 


19 

VU.  HOHNER 


I  he  corresponding  complex  light  distribution  at  the  oulpiil  plane 
I’,  ol  the  processor  lor  each  A  would  he 


g(x.y;A) 


A// 


S<P  P„.9»  H (P„  P.,.11,,) 


<--'p['(P„x  "I  s|„> )  ]c)p„dc|n  .  (5) 

where  (he  integration  is  over  the  spatial  domain  We  assume  that  the 
signal  spectrum  is  spatial  frequency  limited  and  the  bandwidth  ot 
H  tpM.qn)  extended  to  this  limit,  i.e.. 


H(p„.q„).  nt,<'(»<i»,  . 

H<P„-H„)  r" 

0.  otherwise  . 


(hi 


where <»,  =  (Anf  2rr)(p0  f  Ap).und«,  -■ (Anl‘  2rr)(p0  Ap)arelhe 
upper  and  the  lower  spatial  limits  of  ll(pn.t|n).  and  Ap  is  the  band¬ 
width  of  the  input  signal.  I  he  limiting  wavelengths  ol  the  dispersed 
spectra  at  the  upper  and  the  lower  edges  ol  the  filters  are.-'1 


Kf 


Po  +  Ap 

Po-^P 


and 


p,»-^p 
Po  +  -ip 


(7) 


and  the  corresponding  wavelength  spread  over  the  filters  is. 
therefore. 


we  will  describe  is  an  inverse  complex  spatial  filtering  pioeess  We 
consider  the  f  ounei  sped lum  ol  a  blurred  photogiaphic  image  as 

I  (p.q)  Stp.v|)  Dtp.ql  .  (II) 

where  S  (p.q)  is  the  I  oui  lei  sped i  um  ol  unhl lined  image  and  I )( p.q) 
is  the  loiinci  speettum  ol  the  him  tunclion 

In  tie hlu 1 1  mg.  vie  apply  the  him  ied  photogt aplnc  image  It  x .  y )  to 
a  ptesci  ihetl  inveixe  liltei .  as  shown  in  I  ig  2  We  let  the  inv else  Idler 
lunction  be 


hU.y) 

•(*  .y  > 

H(p,q) 

gl»,y);U*,y)--.h(*.y)--sl*,y) 

F  vp.  q) 

_ 

Glp.q^RpflHUp.q)-  Slp.q) 

Fig.  2.  Block  diagram  of  an  image  deblurring  system,  (x.y).  spatial 
domain,  (p.q).  spatial  frequency  domain 


I 

M(p.q)  -  - - . 

D(p.q) 

then  the  output  I  miner  spectrum  is 
(i(p.q)  I  (p.q)  tl(p.q)  -  S(p.q) 


(12) 


(I.D 


AA„ 


4p„Ap 

p.HAp)-’ 


(«) 


If  the  spatial  lrev|ticncy  p,,  ol  the  grating  is  high,  then  the  wavelength, 
spiead  over  the  Idlers  can  he  appio.xmiated. 


which  is  essentially  the  I  ourier  spectrum  ol  ihcunblurrcd  image.  The 
inverse  translorni  is 


g(s.y) 


// 


l(x'.y')  h(x  x’,  y  y  )  tlx’ily '  s(x  ,y ) 


(14) 


4AP 

-*\r  „  V  Po  '  >  ('» 

I’o 

Since  the  complex  spatial  filterings  take  place  in  discrete  I  turner 
spectral  bands  ol  the  light  source,  the  tillered  signals  are  muliiall i 
incoherent.  I  herefore  the  output  light  intensity  tlisti  ihution  is 


v 


Itx.y) 

V 

II  1 

AA(1|g<x.y;An)|-' 

N 

V 

AA„  |s(x.y  ,A11)*h(x.v  .A„)|  . 

(  Mil 

a  t 


where  h(x,y;  An)  is  the  spatial  impulse  response  of  the  tiller  ll(pn.q„). 
and  *  denotes  the  convolution  operation  I  tom  the  above  equation, 
we  see  that  the  white-light  processing  technique  is  indeed  capable  ol 
processing  the  signal  in  complex  amplitude  Since  the  output  intensity 
is  the  sum  ol  the  mutually  meoheienl  nation -band  uiadiaiiccs,  the 
■  annoying  cohcicnl  ail  duel  noise  can  be  suppressed 

In  addition,  the  white-light  source  contains  all  the  visible  wave¬ 
lengths  ol  I  fix-  elccf  lotuagtictic  wave,  ami  theteloM-  is  p.nli,  nlailv 
suitable  lor  color  image  processing.  We  luitlier  note  t hat  the  while- 
light  processor  of  l  ig.  I  can  also  be  used  lot  coherent  and  partial 
cohcicltt  light 

V  PHOTOGRAPHIC  IMAGK  DKRI.t  WRING 

Oncol  the  interesting  applications  ol  optical  inlormation  processing 
is  restoration  of  blurred  photographic  images.  I  he  dchlin  iing  that 


which  is  the  imhlurrcd  image  s(\,y  i 

We  note  that  photogiaphic  image  dehlurimg  by  coherent  optical 
pioccssmg  technique  was  illustrated  by  I  sii|iiicln-”  m  I9(i7  I  he 
inverse  spatial  liltei  was  sy  nthesi/ed  hy  the  combination  ol  an  ampli¬ 
tude  and  a  pure  phase  filter  I  his  can  also  be  accomplished  by  a 
holographic  synthesis  technique.  I  he  preparation  of  such  a  phase 
filter  bv  hologiaphic  technique  had  been  studied  by  Stroke  and 
/ech,  '*  ami  by  I  ohm. inn  and  I’ai  is.-"1  Ncvei tireless,  the  llologi apltic 
synthesis  (eclmiqire  also  sullers  one  disadvantage,  namely  a  low 
dilliaction  vlluienvv  Meniioit  must  also  be  maile  ol  ilic  image 
ilebhiiiingwillicompiilei  generated  phase  liltei  obtained  In  I  sii|iiichi. 
Honda,  and  I  ukava  I  hw  had  show  nv;n  ions  ellects  due  to  .impli- 
luvlv.  pli.isc.  and  aiuphliivle-phase  lilleiings  Auothei  ililciestmg 
result  obtained  Ivy  Horner"-’-’  should  also  be  mcntiimed  lie  has 
shown  lh.it  optimum  image  ilchlurring  may  be  obtained  with  least 
mean-square-eiioi  I  liter  mg. 

We  will  bi icily  describe  the  synthesis  ol  a  phase  liltei  that,  when 
combined  w  ith  an  amplitude-  filter,  can  be  used  lor  image dehlurring. 
I  oi  es.imple.  Ivl  the  amplitude  liuiiMimtuiH.-v-  ol  a  liiic.n  smeaievl 
point  image  be 


Kv) 


I  I 

I.  ,  Ay  y  ,  Ay  . 

I).  ol hei wise  . 


(ISI 


where  Av  is  the  smcai  length  II  the  him  ml  image  is  mseiled  in  the 
input  plane  I’,  ol  a  coherent  optical  •  loccssor  ol  l  ig  I .  the  complex 
light  livid  on  the  spatial  licvpieney  1  .me  would  he 
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stnlqAv  2) 

Hi))  -  Av . .  (If>) 

qAy  2 

(he  familiar  urn.’  I  unci  ion  In  p  nimplc.  I  Ik  him  nil  image  may  be 
dcbluiied  hy  invcise  1 1 1 1  c  ■  i  n  ^  I  mlu  mi  we  lettbe  litlei  he  a  coinbma- 
lion  ot  an  amplitude  and  a  phase  Idler,  as  shown  m  l  ip  I  he  filler 


Fig  3  Inverse  spatial  filter  for  deblurring  a  linearly  smeared  image  (a) 
Amplitude  filter  function,  and  |b)  phase  filter  function 


<i(q)  ,  I  (e|>  AM  I  4  -  fill))  H(e|)  e\pll>(,  1)1 

i  I  tql  H*(ej»  exp(  iy„ql]  (20) 

It  is  clear  l  hat  i  lie  In  si  lei  in  is  the  spectrum  due  to  the  amplitude  filler 
alone,  which  w  ill  he  dill  laded  on  I  he  optical  axis  at  the  out  pul  plane 
I’,  ol  I  tg  I  the  second  and  thud  leims  are  the  lesioied  I  miner 
specli a.  in  which  the  lesioied  images  will  he  respectively  ditlracled 
around  >  >(1  and  y  \  al  the  output  plane  I’, 

An  ink-testing  alternative  approach  related  to  holographic  syn- 
thesis  ol  the  Idler  has  recently  been  reported  h\  Vasu  and  Kogeis  14 
We  now  briefly  describe  a  while-lighl -processing  technique  lor 
smeared-iinage  debliu  ring  We  place  a  smeared-image  transparency 
in  contact  w  nil  a  sinusoidal  phase  grating  at  the  input  plane  l’|  ol  the 
white-light  pi  ocessorol  fig  I.  I  he  mathematics  ol  this  base  already 
hecn  described  in  Sec  2  for  simplicity.  we  assume  that  the  input 
transparency  is  spatial  frequency  limited  and  that  the  smearing  is  in 
the  y  direction  ol  the  input  spatial  plane.  We  further  assume  that  a 
narrow-hand  dehlurnng  filter,  for  a  given  wavelength  An  as 
described  ealliei.  is  provided. 

II  we  insert  the  dehlurnng  filter  of  ll(/f)  over  a  narrow  spectral 
band  ol  the  snieared-signal  spectra  l  [n  (At  2:r)pi,./l]at  A  -  A0. 
then  the  complex  light  field  at  the  output  plane  1’ lor  a  given 
wavelength  A  over  the  dehlurnng  tiller,  can  he  evaluated  b>  the 
following  inteeiul  ei|iiation 
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glx.y .At 

(  r 

1  J  J  "* 

"t 

2  n 

exp  1  (i x  t  /Iv) 

At 

cl  < *li  fi  . 

(21) 

where  C  is  a  complex  constant  and  II (/it  is  the  dehlurnng  filter  The 
corresponding  output-light  intensity  distribution  can  he  approxi¬ 
mated  bv 


l(x.y)  J  Igtx.y  ;A)l'd A 

AA 


k  AA|l(.x.y)  exp(ip0x)  *  h(yl|‘  . 


(22) 


transfer  function  is 

ll(q)  A(q)  e\p|n/,(q)  | 


(l?l 


We  let  I  he  the  minimum  tiaiisiiiiltance ol  the  amphiiide  lilici 
and  the  relative  degree  ol  restoration"’  “  is  delined  as 


l> 


I _ f  f  (q)H(q) 

l„,Aq  J  Ay 


dq  V  KKI'  ,  , 


OKI 


where  Aq  is  the  spatial  haiulw  idtli  of  the  I  diet  We  note  that  pci  led 
test  mat  ion  is  mint  tamable,  since  it  icqimrs  that  I  appio.n  lies  /cm 
*  I  hete  is  also  anothci  shoilcoumig  loi  lowei  I  it  would  result  in  a 
lowei  signal  to  noise  nitio  llv  taking  the  an  omit  ol  him  gi. mi  noise, 
.i  ill  ill  mill  ill  value  ol  1 1|(  with  a  lau  Iv  good  signal  In  noise  lalio  ill. IV 
.  he  attained 

We  now  assume  Dial  the  phase  Idler  is  a  holographic  iv  pc;  that  is. 


Ilq)  ,  I1  4  eosfcMq)  t  y„q|  | 


(I'll 


where  v„  is  an  arbitrary  constant.  I  hen  the  tillered  I  miner  speed  um  is 


where  AA  A^tdAp  p,,)  is  the  narrow  spectra!  band  ol  light  over 
the  dchhiiiing.  Iiltet,  h(\  t  is  the  spatial  impulse  tesponse  ol  lll/i).* 
denotes  the  convolution  opeiation.  and  k  is  a  piopoitionahlv  con¬ 
stant  I  tom  the  above  equation,  we  see  that  the  liluicd  image  is 
dilliaclcd  atouml  the  optical  avis  al  the  output  image  plane  I’, 

As  an  cvpciuuciit.il  illustration.  I  ig  4(a).  shows  the  liricar- 


41a) 


OPTICS  OPTICS 

4(b)  4(c) 


Fig  4  Linearly  smeared  image  deblurring  (a)  Smeared  ob|ect.  (b)  da- 
blurred  with  whito  light  technique,  and  (c|  deblurred  with  coherent  light 
(Zhuang.  Chao,  Vu,  1981”’) 
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Fig  6  Smuarml  image  diibliirrmg  m  color  In)  Bloriml  imugu  and  (b)  do 
blurred  with  white  light  technique  (Yu.  Zhuang.  and  Chao,  to  be 
published1*) 


sr.icjred  photogidi'liii.  itn.igc  ul  ,i  «mil  n/>in  i  .iv  .1  bliiim)  i>b|wt 
l  igurc  41  h  I  shots  \  I  he  tic  bin  tied  im.igc  nbl.uiu-d  vv it  It  the  w  lute  -  light - 
piviwwiig  kvlmk|tic.  .md  I  ig  4(c)  slums  thv  usiili  . il'i.mud  with 
the  cohciem-|irou-"int!  technique  I  \pheit  details  .tie  given  m  Ket 
.1?  I  rom  these  lestilts,  we  see  that  lhe.11td.1cl  nuise  w.is  siibsi.uiti.ilh 
reduced  with  the  whilc-light-processing  tecbmque  I  he  debluned 
image  obtained  with  the  eohetettl  leehnu|iie  appeals  10  be  sh.it pet 
than  the  one  obtained  with  the  white-light  leclum|iie  because  ol  t lie 
high  spatial  eohetenee  ol  the  light  souiee  llmcevei.  the  <li.iwb.uA 
can  be  osercome  it  one  uses  a  smallei  white-ligbt  somce  lie  .  a 
pinholel  and  a  broad-speetial-band  deblumng  Idlei  to  covet  the 
entne  snieaied  signal  spevtiuiu  lUsipn  ul  mu  h  a  bm.nl  band  tail 
shaped  liltci  is  tindei  investigation 

Since  the  white-light  sou  ice  contains  all  v  isible  wavelengths  ol  the 
specti  mu.  the  w  lute -fight  pioccss  is  suitable  loi  coloi  image  deblm- 
nng  "■  I  et  us  place  a  smeared  color  phologtapliic  image  at  the  input 
plane  I*,  ol  the  white-light  pioccssoi  ol  lag  I  At  the  spatial  Itetpieti- 
c>  plane  I’,,  two  sets  ol  I  ourier  spectia  are  snieaied  in  lainbow 
colors  In  color  image  deblurring,  we  place  three  prunaiv  color  (1  e  . 
red.  green,  and  blue)  narrow-spectral-hand  dcblumng  fillets  m  the 
apptopriate  locations  in  the  spatial  lievpiencv  plane  Since  the  de 
binned  images  hum  each  ol  (lie  piiuiaiv  coloi  Itlteis  aie  unit u.ills 
incoherent,  the  output  image  niadiance  will  be 

Its.y)  -iA  f  |g(x  .v  .  At  )|"  I  |  g(  s .  v ;  A  g  1 1 

I  |g(.\.y;Ahl|-’  ]  .  (2b 

where  AA  is  the  narrow  spectral  band  of  the  debfurrmg  filter; 
g(x  ,y  ;A)  in  the  dcblurrcd  image,  and  Af.  A^,  ami  Ah  ate  the  led.  green, 
and  blue  primary  color  wavelengths. 

I*  or  an  experimental  tihisiiation.  a  linear -smeaicd  color  photo- 

. . .  •«!  1  M.im  1  inm,  /  . . .  n 


Fig  6  (mage  subtraction  with  encoded  extended  incoherent  source 
MS.  multiple  slit;  P2.  input  plane,  P3.  filter  plane.  P4.  output  plane.  Ul, 
12,  L3,  convergent  achromatic  lenses  of  focal  length  F  (Wu,  Yu.  to  be 
published44) 


Fig.  7  Image  subtraction  (a)  Input  objects,  (b)  output  with  incoherent 
technique,  and  (c)  output  with  coherent  technique  (Wu.  Yu,  to  be 
published  44) 


graphic  image  ol  a  word  “optics”  is  show  n  m  I  ig  5(a)  I  he  dehlurred 
color  image  obtained  with  the  white-light  technique  is  shown  in  f  ig. 
5(h).  Strictly  speaking,  a  Ian-shaped  dehlumng  lilter  o\ei  the  entire 
smeared  color  !  oimei  spectra  should  Ik-  used  Such  a  lifter  is  under 
investigation 

4.  IMIOHH.RAIMIK  l\1  \(.Y  SCRTH  ACTION 

Anolhei  micicsimg  applieation  ol  optical  photogiaphic  image  pro¬ 
cessing  is  image  subtraction  image  suhliaction  may  he  ol  value  in 
main  applications  such  as  ui ban  development,  highway  planning, 
earth  resources  studies,  remote  sensing,  meteorology .  automatic  sur¬ 
veillance.  inspect  ion  etc.  Optica!  image  subtraction  may  also  he 
applied  to  tommumc.Uious  As  a  means  ol  h.mdw  i« h !i  compiessroit, 
loi  example,  it  would  he  necessary  to  liansmil  onls  the  dillciences 
between  images  m  successive  cycles,  ralher  Ilian  the  entire  unage  in 
cai h  cvelc 

Optical  image  synthesis  h\  complex  amplitude  subtraction  was 
desciibcd  by  (i.thor  el  a  I  1  I  he  technique  involves  successive  iccoul- 
iugs  ol  two  or  more  lompfc.x  dilfiaction  patterns  on  a  hologiapluc 
plate  and  the  siihsequcni  teprodueuon  ol  the  composite  hologram 
images  A  few  years  fatcr.  Mromfey  eta)  u  ‘“'described  a  hologr.rphrc 
I  ouiiei  subtraction  technique,  loi  which  a  real-time  image  and  a 
prevumslv  recorded  hologram  image  can  he  subtracted  Although 
gooil  image  subtraction  by  then  expetimenis  was  reported,  it 
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I'piu.il  pioeessot  i*l  I  10  I  \l  lliospat. a!  Iroijisoin'%  piano  P  .  tho 
complex  l.ehf  iliNtnhui ii»n  Io:  oaoh  waxelcneth  a  .m  tho  heht  somcc 
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f  iq  8  Ituf’h*  .%  Imttmly  ipm  ml  in  tog  intensity  (a)  Originul  imogo.  (h) 
halftone  iiorty  df  ul  (i.)  isophot  image  iStraml.  1  li  7b’  ) 


lrc\piciK\  plant  Mnu’  iho  dispersed  I  ouner  spootia  varies  with 
ii  >1'  it  I"  tht  \\  a  \  e  le  i  ct  h  of  tiit  hohi  mmiiio  wr  wmilil  noon  a 
tali  s  I  la  pi  1 1  in  at  I* »  v  *  •  npi  loate  tho  w,i\  ■  leiipl  Ii  \  anal  ii<ii  Wo  lot 
t  ho  lafi  shap»  il  oi  at  ilu*  ho 

Mh|l  |  I  sui/<t,i||  i.Wj 

ami  tho  mil  put  .lilhai  I  i»*n  «  an  hi  nImiw  ii  a 

0 1  X .  \  \  ^n|n(|\.1  /<tt>  s,t\.N  •  /<t| 

*  I  I  2)  (  'N  !  |\  .  V  I  N.l  X.V  I 

S,(\.\  ♦  s , ( \  \  •  .'/*„)  |  I 

o\pnp,,\i  iJ?h) 

I  hits,  tho  sii  ul  r  at  I  iimi  "I  tin  iwh  input  -.ipnal  .  |  •.  (  i  \  \) 
s  \  .  v  )  |o\pl  i  p>(  \  )  tail  ho  st  on  at  iho  opt  it  a  I  axis  ot  t  ho  •  ml  pul  piano 
\Vo  point  out  that,  in  practice.  it  is  if  1 1 1  u  nit  to  obtain  a  w  Into  ho  lit 
point  source  Howovci.  tins  shiMUoimn;  o.tti  ho  ouimnu'  with  a 
soui oo  encoding  technique  that  \\e  will  dcMirhe  “ 

Smoo  nn.ipo  siihtiatlmn  is  a  mu  limit  in  tonal  (I  Ml  pintovsinp 
operation.  msto.nl  ol  a  pinnt  Minito  o|  held.  a  lino  son  roe  ol  heht 
pci  pcndu  ular  !o  tho  directum  ol  the  »wt»  in  rapes  oan  ho  utilized  In 
addition,  tho  collet  cnee  lequuefncnl  lot  mi.tpc  suhti.ntion  in  o.iiK  lot 


Fig  9  Pseudocolor  encoding  with  halftone  screen  (a)  Obtained  with 
cohimmt  souich  and  tt>)  ol>tmnmJ  witii  white  light  source  (Tm,  Yu  and 
Chen.  19/04) 


c\ci\  point  pan  holwccn  the  two  nnapes.  and  a  stnot  tolioictUC 
i Oi j i in i’ii io t it  is  not  in  nlnl  In  oi  hoi  w oi ds.  it  is  possihlo  to  oik  ode  an 
i  x  lom  h  1 1  i  iiv  <  1 1 1>  1 1  ul  m  mu  o  1 1  >  ai  hi*  \  i  Hi  poi  ill  pa  it  i  ol  it  i ,  me  io 
t| m lenient  I  o  do  so  wo  appls  the  V  an  C  ittoi t  /omiko’s  ihoorom  J‘ 
\  point  pan  eoheicnl  lunclion  can  ho  deter  turned  \n  a  result  one 
can  ciuodc  .in  cMcnded  meohetenl  soiiuo  with  a  n  m  it  is  1  it  mask 
W  it  h  i otri  oi H  o  t .*  I  io  f*.  the  i  ohot  cut  lum  t  nm  at  llu  itip  u  plain  P  .  is 
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where  d  is  ihc  spacing  nt  the  slits,  s  is  the  slit  vs  idth.  and  N  is  the  total 
ihi m her  ol  slits  ovei  the  extended  sonice  I  he  degice  ol  spatial 
coherence  tor  each  pair  ol  points  sep.ualed  hv  a  distance  'h  t  is 
t  tier  el  ore 


m(2H  t  Si: 


l  rom  the  aho\e  equation,  it 


a  high  degiee  ol  point-pan  coheieiice  van  he  achieved  In  coinpai  ison 
with  a  single  slit  case,  the  mtensitv  ol  the  output  image  can  he 
increased  N  told  In  olhei  wotds.  the  sou  ice  encoding  pi  ovules  an 
etlicient  uttli/ation  ol  the  light  source 

In  our  expenmenta!  demonstration,  a  nmltislit  mask  with  a  25  ^m 
spacing  and  2  5  width  was  adopted  tor  the  souice  encoding  I  he 
tight  source  was  a  mercutv  arc  lamp  with  a  green  tiltei  I  iguie  "'la) 
shows  two  image  tianspatencies  as  input  ohiects  I  iguie  ‘Ah)  shows 
the  subtracted  image  obtained  with  the  incohcicnt  technk|ue  In 
comparison,  we  also  prov  »de  the  subtracted  image  obtained  w  it  h  the 
coherent  technique  as  shown  m  l  ig  “*tci  I  tom  the  results,  we  see 
that  the  one  obtained  with  the  incohcicnt  technique  oilers  hettei 
image  quuhtv  and  less  attilact  muse 

I  mails .  we  w otild  stiess  that  the  source  encoding  ma\  be  extend¬ 
ed  to  white  tight  processing  technique  and  a  ptogi.im  is  cmientlx 
under  wav  to  mvestieate  this 


giaptiic  dv  n.mik  i  an  ge  v  o:u  pi  ess  ion  pi .  tpu'cd  hv  Mue  tie  i  arid  (  aul- 
tield  '  I  tu  I iavc  shown  that  mditi.u  v  photogiaphu  emulsion-  can 
lecoid  ml. um.it u *n  ol  high  dvu.imu  range  that  we  usu.illv  do  mit 
lecogm/c  1  la  \  u  .i  d  a  ■  i.rl.ur  rclul  U\  liimpie  that  Ua-  little  ellvci  on 
the  d  *.  an«l  love  -patuil  hequeticV  l  out  cut-  *-!  t tic  p t ; >  M. .gs aphic 
image  I  he  icliet  patiem  is  more  like  a  -patia:  gradient  o!  the 
i  c\  or  ilcd  d<.  ns  it  \  pall*.  i  ti  -«u  h  I  iia!  the  -t '  i  la.  c  h.»s  t  he  c: .  .\\  -t  i  el.et 
w  heie  t  !ie  -pat  ia I  dciiMt  v  patten,  c  hang.-  r:.« *-r  i  .ip  d  .  \ .’  t  .  -ugl.  r  tie 
relict  in  .age-  were  i:.*t  readilv  observed  wi'.r.  dith.-e  i  .  iiiat  i.  >u. 
using  a  m».t!  ii.votierent  pto^es-mc  ‘•t.u.ae  the.  were  .P»!c  to 
sliovv  a  w  ide  d v  iiamu  lance  .*!  (Pc  pimtog:  aphu  -mage  a-  deputed 
m  f  :g  l<)  I  he  tec  fi usque  was  applied  r*.  td.u  ►  mg  the  d a  ».< mi. portent 
ill  I  lie  I  i  m:i  ic«  plane  o I  the  lei  i  of  clU  * '  .ul  h/l.'  I  *  •  nti  t  lie  le  he!  of  lace 
ol  t lie  pfnMogi .ipt' a  titm 

6.  SI*A(  l-VAKIVM  I*RO(  KSSIM.OI  1*110 KK.KAPIIK 
IM  \(.l  S 

In  t tie  prev  .ous  sections  we  have  disc  u--cd  on!>  the  space-invariant 
processing  operation  ol  photographic  image-  In  other  word-,  each 
image  point  to  be  processed  i-  allccled  bv  the  same  pr.ue-- g 
operation  Although  the  -pace-v atiatil  processing  concept  ha-  been 
Known  ami  succc--lu!lv  applied  in  -v-ten.  thcoiv  ami  digital  process¬ 
ing  ted  i  tuques  '  the  application  ot  th.s  w  i  Mitcptlo  optical  processing 


5.  NONI.1NKA K  l*llO(  I  SM\(.  Ol  IMIOKM.K  \IMIH 
IMACIS 

So  lai  ill  pi  ev  tons  sect  lot  is  w«liav  c  »!i-c  u-  cd  Imc.ii  -  pal  .al  -  mv  ,u  i.ml 
operations  I  liete  are.  however,  techniques  available  tot  nonlinear 
processing  operations  u  *  One  such  appioadi  t\  an  opinal  Immo- 
motphic  lilteung  svstem  which  has  been  quite  simc-hith  applied 
Using  dig  it  til  techniques  '*  l<  c»  c  n  1 1  v  .  Kato  a  in  I  (  on  idm.m  •  ju  ujiu  .cd 
nonlineai  processing  ol  photographic  images  with  a  ct./rent  optical 
svstem  using  h.tlllone  scieeit  ptoecsses  I  he  extension  «•!  kato  and 
( i  ood  mail's  ha  It  tone  sci  ccn  U\  Imiqiic  to  iioiinionotouic  .11  id  mono 
tonic  nonlinear  image  piocessing  vv.is  stihscc|uenii\  pm -ucd  bv  Saw 
chwk  and  I  )»islueIl.Vl  Si  1  and.* '  as  show  n  m  I  ig  .v  and  I  m  (mod 
man.  and  (  dun  *■  Application  ol  the  halltone  scicen  U\  liniques  loi 
pseudocolor  encoding  ol  monotone  images  with  coheirnt  opiic.il 
processing  was  also  icpoiied  bv  I  111  and  (modm.iii'*  and  with  an 
incoherent  source  hv  I  ai.  N  11.  .ind  (  lien  **  as  shown  in  I  ig  d 

We  note  one  nia|ui  limitation  with  the  lulltoce  scic.nmg  tech- 
nique.  the  spatial  resolution  ol  the  photographic  image  i-  limited  hv 
the  halltone  sv  teen  iiselt,  w hu  ll  is  gcncialiv  sevcial  01  «h  is  ol  magni¬ 
tude  lowei  than  the  pliotogi apluc  him  It  might  be  possible  to 
piodiice  .1  halltone  pidiiie  with  icsoluliou  up  to  |0(i  I  mm  using  .1 
computer -ill  iven  wnlei  but  the  pi  mi  would  he  lomplcx  and 
cosilv 

1  he  nonlineai  processing  technique  can  also  be  .u  lncvcd  with  a 
I  abiv-l’ciot  inter  lei  omet  lie  metliod  .is  icpoitcd  bv  li.n  tlioiomcw 
and  I  ee,'m  and  with  the  utilization  ol  him  nonlinear  Hv  icpoiied  h\ 
lai,  C  heng,  and  >u  *'  Mciilmn  must  be  made  that  leal-time  non¬ 
linear  processing  ol  photogiaphu  images  v\  it h  liquid  crvstal  v  alv  es  is 
cuiientlv  being,  puisued  bv  Mnhaelsou  and  S.iwdmk.'  and  bv 
Sol  lei  el  al 

It  is  possible  lo  use  the  minimi  him  iionlincai il v  in  peilmm 
logarithmic  ami  exponential  lianslm malmiis  t’sing  the  basic  II  I) 
curves  ol  various  hints,  carelullv  choosing  the  appiopnaie'y  and 
using  two  step  v  out. 11 1  pi  lilting,  Vu  .11  lueved  a  logai  it  In  tin  amplitude 
transmittance  loi  an  input  exposute  i.inge  ol  Ml  I  lu  inverse  Hans 
tormaliou  ol  the  hltcrcd  logarithmic  signal  can  also  hi  pethiimed 
using  a  two-step  approach  I  he  details  and  specific  lilms  used  can  be 
found  in  Rel  59 

l  inallv.  we  note  that  theie  is  an  iiileiesimg  tcchnnpic  ol  photo 
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Fig  10  Dynamic  range  expansion,  (a)  Overexposed  photograph  and  (bt 
bright  held  nitririvnl  from  iho  rnliol  image  of  the  original  negative  (Moeller 
and  Cmillmld.  19B0hl*) 
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in  relatively  icccni  In  I OnS,  t'lilion.i'-  proposed  .1  gciici.il  optical 
spacc-vaiiant  processing  concept.  analogous  to  the  I  miner 
transtoim  technique,  in  lei  ms  ot  ssstem  eigenfunctions  I  he  leeli- 
mque  is.  howcict.  nl  onl>  lheorclic.il  interest,  since  n  is  not  l  noun 
how  to  Inul  the  set  ol  eigenfunctions  Attempts  had  been  made  hs 
Hem/.  Allman.  and  I  ee  in  optical  spacc-s.ii i.mt  piiKcssmg.  using 
Hiatus  inn  In  plication  Some  calculated  icsulls  weic  1  c  ported  in  t  hen 
paper  In  .1  recent  article.  Deen.  Walkup.  and  llaglci'"1  reported  a 
sp.ii.c-s. 111. ml  opei.ilion  using  volume  linlogi.ims  I  lie  most  gi.iplnc 
example  ol  space  s  .11  lain  pi  ocessing  is  the  geomcti  ic  tiansloi  mat  ion 
b\  Hi s ngdahl.' ‘  who  useil  a  computer-generated  hologram  lor  the 
operation  Mention  must  also  be  made  ol  Sawchuk's  work.'"  using 
the  computet  technique  ot  space-sanant  image  restoiation  bs  coor¬ 
dinate  traiislot  illation 

We  will  now  Jesetibe  a  basic  concept  to  aclnese  optical  space- 
sanant  processing  similar  to  that  ol  homomophic  piocessing  illus¬ 
trated  in  the  pres  ious  section  A  genet al  space-sanant  system  mas  be 
desenbed  bs  a  block  diagram,  as  shown  in  I  ig  Ilia)  I  he  mput- 


Fig  1 1  Block  diagram  of  space-variant  processor,  showing  nomenclature 


output  relation  mas  be  described 

g(x.y)  ^^hfx.y  .xV)  Hx'.y'fdx'ds'  .  ( 

where  h(x.y;\'y ')  is  the  space-sanant  impulse  response  Needless  to 
sas .  il  hls.s.sVI  is  described  as  a  I  miction  til  the  coot  dilute  dillci- 
encc.  ic.  h<\  x'.s  s'),  then  the  system  is  space-insariant  (or 
example,  a  technique  ol  optical  space-sanant  processing  imp  be 
accomplished  bs  the  block  diagram  of  I  ig  ll(bl  In  other  words, 
w  hones cr  a  space-variant  process  can  he  decomposed  .is  in  I  ig 
I  1(b),  it  mas  be  performed  by  the  necess.r  s  geometric  translorma- 
tions  I  lie  geometric  transform  signal  can  be  pioccsscd  hs  a  space- 
insariant  ptocessoi  and  then  the  required  mseised  geoir.cinc.il  tians- 
I01m.it  ion  pci  loi  tiled 

In  a  moie  iccciil  aitnle.  Casascnl  ami  I’salns'-'  base  shossu  that 
geometnc  tianslormalions  can  also  be  accomplished  bs  means  ol  .1 
nonlinc.il  beam  st  .inning  ties  ice  which  wi  lies  the  il.iiu  miln  .111  .>pli 
cal  light  salse  Although  this  method  sacntices  the  pa  1 .1 1  Ic- 1  piocess- 
ing  capability,  it  can  process  in  near  leal-time.  Although  mleiest  in 
space-sariant  processing  ot  photographic  images  is  hurls  recent,  in 
slew  of  this  interest,  lunhcr  progiess  can  he  expected  in  future 
research. 

7.  AKdllVAl  SIOK  \(;K  ANI)<  Ol  OKl  MIANf  I  MI  Nl 

Arclusal  slot  age  ol  color  lilms  has  long  been  an  mesolscd  pioblem 
bn  him  11  id  11st  1  ic.  I  lie  1 1 1.1 1 01  1  cast  111  is  1 1 1,1 1  I  In-  01  gun  it  lives  used  11 1 
the  color  hints  ale  iisiiails  mist  able  iintlci  pioloiigcd  si  01  age.  and  I  Ills 
causes  a  gradual  color  lading  Although  there  an-  seseial  as.nlable 
tccliniqucs  lot  pieseising  the  color  images,  all  ol  them  possess  cci 
lain  drltmlc  draw  backs  With  one  ol  the  most  cominonls  used  tech¬ 
niques.  that  ol  icpclitisc  apphc.it  ion  ol  |ii  imai  s  coloi  I  liters,  one  cvin 
preserve  the  color  images  in  three  separate  rolls  ol  black  and  while 
him  lo  rcpioduce  the  color  image,  a  ssstem  with  three  primary 


coloi  pioicctois  must  be  used  I  hese  lilms  should  he  projected  in 
pel  led  unison  so  that  the  piimars  color  images  will  be  precisely 
iccorded  on  a  Itesh  10II  ol  color  him  However,  this  technique  has 
two  maim  disadvantages  (I)  the  storage  volume  lor  each  him  is 
tripled  and  (2l  the  icpioduction  ssstem  is  rather  elaborate  and 
expensive 

In  this  section,  we  will  describe  a  white- light  processing  technique 
loi  aichisal  storage  ol  color  hints  I  his  technique  may  be  the  most 
elhcrent  technique  exisimg  to  dale  I  his  technique  also  allows  direct 
s  lew  mg  capability .  and  may  he  suitable  loi  library  applications 

I  he  use  ol  monochrome  transparencies  to  retrieve  color  images 
was  Inst  report ed  by  Ives  in  1906.  He  introduced  a  slide  viewer  that 
produced  color  images  by  a  diltracnon  phenomenon  (iratings  of 
either  dillcicnt  spatial  frequencies  01  u/imuthal  orientation  were 
used  More  recently .  Mueller'"' described  a  similar  technique  employ¬ 
ing  a  incoloi  grid  screen  lor  linage  encoding  In  decoding,  he  Used 
thiee  qiiasi-monochromatic  sources  lor  color  image  retnesal  Since 
then,  similar  woik  on  color  image  retrieval  has  been  reported  by 
Macovski.  "  (ooiisson  and  Kinaris.  1  and  Yu 

We  describe  the  technique'''  -  loi  spatially  encoding  the  color 
in  1  1  illation  Horn  a  single  color  image  onto  a  single  black  and  w hue 
(HW  )  liume  I  Ins  is  done  by  putting  a  color  filter,  a  Ronchi  grating, 
the  color  frame,  and  the  unexposed  HW  film  together,  and  exposing. 
A  tuple  exposure  is  made  onto  the  HW  film,  changing  the  color  filter 
(red.  green,  and  hlucl  and  the  angular  position  ol  the  Ronchi  grating 
each  time 

I  hi  example,  the  first  exposure  is  made  through  a  red  filter  with 
the  giat mg  at  0  I  he  second  exposuie  is  made  through  a  blue  filter 
with  the  grating  oriented  by  60  ,  and  the  last  exposure  lakes  place 
through  a  green  filter  w  ith  the  grating  at  I  20"  II  the  three  recordings 
are  pinpeily  recorded  on  the  photographic  film,  then  we  have  a 
multiplexed,  spatially  encoded,  black  and  white  transparency. 

II  we  place  the  spatially  encoded  transparency  at  the  input  plane 
I* j  of  .1  white-light  processor  shown  m  lag  I.  but  without  the 
diffraction  grating  ( (XI.  the  different  orders  of  the  image  .spectra  are 
lineal  Is  dispersed  111 1  a  in  boss  colors  with  respect  to  the  *.  and 
axes  II  the  spatial  l  requeues  ol  the  grating  p0  is  assumed  sufficiently 
high,  then  the  orders  of  the  smeared  color  image  are  physically 
separated 

I  o  letnese  the  color  image,  we  allow  three,  first-order,  smeared 
spectia  to  pass  respectively  through  a  red.  a  blue,  and  a  green  color 
filter,  as  shown  in  lig  12.  Uteri  three  mutually  incoherent  primary 
coloi  images  will  be  super umpo  ed  al  output  plane  I  bus,  we  see  that 
a  multieoloi  image  is  reproduced  hs  the  white-light  processing 
technique 

Figure  I A( a >  shows  a  retrieved  color  image  obtained  with  this 
technique  In  this  experiment,  we  used  only  red  and  green  color 
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Fig  12.  Color  image  retrieval  spatial  filter  (Yu.  1980”) 
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Fig  13  Color  <mayo  rcrccmlod  by  multiple*  tochrmjur>  on  I3W  film  (a) 
Original  image  and  |b)  retrieved  image  (Yu.  1980") 


liltcrs  tor  i he  spatial  encodings  I or  comparison.  \u'  provide  the 
original  color  picture  as  shown  in  I  rg  I  hbf  We  see  that  ihe  color 
reproduction  bv  this  archival  stoiagc  technique  is  mipiessivclv  l.uih 
lul  and  i he  resolution  is  good. 

An  eMension  ot  the  above  whitc-huhi  technique  can  be  used  to 
restore  or  enhance  ihe  colors  m  tailed  coloi  tilms  '  l  sing  a  standard 
coherent  optical  processor .  Hornet  ‘demonstrated  a  similar  solution 
tor  the  Apollo  coloi  hints  degraded  1>\  i.nli.ition  lit  oulei  span 

8.  PSKI  IXK  Ol OK  I  M  OOIM.  OI*  IMIO  TOOK  AI’IIU 

imaois 

Most  ol  the  optical  images  obtained  m  \  ai  ions  seient  ilu  applications 
arc  densttv -modulated  black  and  white  images  Human  ohscivcrs. 
however,  can  perceive  van.itions  in  odors  hi  tter  than  gras  levels 
I  bus  a  color  encoded  image  can  olteti  pmv  ule  heller  v  istra!  ihsci  imma- 
fion  INeudocolor  em  oiling  b\  computer  lei  liuique  has  been  widclv 
used  m  applications  whcic  the  images  aie  imtiallv  digitized  While 
the  computet  technupie  is  the  logical  choice  toi  digital  images, 
opt n  ,tl  ptot  essiug  lei  hiu<(ues  mav  he  nioie  .nlv.inlag*  oil .  Joi  apph 
cations  where  the  initial  images  aie  analog  photographie  images  mu  h 
as  aerial  or  vtav  photographs 

Histoi  icallv  the  use  ol  pseudocoloi  encoded  lilteis  m  an  imaging 
system  was  lust  mtioduced  by  Khcinberg  in  IX9f»  lie  icpotled  the 
apptre atmii  of  color  Idlerrng  in  rnuroveopy  He  showed  lhal  bv 
simple  color  spatial  tillering  he  was  able  to  enhance  small  details 
Hutch  '  applied  color  spatial  tillcis  lor  enhancement  ot  m  (dilation- 


Fig  14  Pseudocolor  encoded  x  ray  image  (Chao,  Zhuang.  Yu.  1  98080) 


dependent  and  high  spatial  licqucucv  details  In  a  more  leeent  arti- 
i  le.  Hes».os  .i ml  Shand  pioposed  a  coloi  encoding  technique  using 
an  extended  polvchi oin.il .c  light  souiee  I  lie  proposed  technique 
encodes  the  image  hv  spatial  licqucucv  lailiei  than  hv  duisilv  N  u  el 
a  I  have  mttodiu  ed  an  alternative  appto.n  ii  using  the  one -step  i  .im¬ 
ho  w  hologiuphic  eiuodmg  technique,  w  heie  the  encoding  was  also 
done  m  spatial  frequency.  Ihe  encoded  holographic  image  can  be 
reproduced  bv  simple  white-light  illumination  Irulcbetniiw  '*  ?mro- 
duced  a  w  lute-light  pioeessmg  technique  encoding  the  image  density, 
in  which  the  pseudocolor  images  weie  obtained  thiough  color  t titer¬ 
ing  at  the  spatial  licqucucv  pi. me.  I  in  and  (lOodmaiT'  dcscnhcd  *l 
technique  using  a  special!)  luhi  icated  halllonc  set een  to  obtain  den- 
sitv  pseudocolor  images  with  cohcrenr  optical  processing  lar.  Yu. 
and  (  In’ii  1  adupied  then  halllonc  liatniiqiie  to  geneiate  pseinlo- 
color  densitv  images  with  a  w  lute-light  soincc  However,  with  the 
hall  lone  screen,  a  immhci  ol  disci  etc  lines  due  to  image  sampling  arc 
genciulh  present  in  the  color-coded  image  1  has.  small  details  and 
low  contrast  features  ot  the  image  can  be  lost  m  the  haltlone  tech¬ 
nique  A  new  densitv -pseudocoloi  mg  technupie  thiough  contiast 
reversal  was  reported  recenllv  In  Santam.ma.  (ica.  and  Hisois 
Although  this  technupie  oilers  seveial  advantages  over  pievious 
techniques,  the  optical  sv stem  is vpute  elaborate,  and  it  icquitcs  both 
nit  oiu  i  nil  and  •  oheo  ul  stun  ces  I  or  ihe  j  iscudm  oloi  mg  I  let  a  use  a 
eoheienl  sotme  is  used,  the  eoheient  ailitaet  noise  mav  not  be 
avoided  Keci  ntlv.  A  u  desciihed  a  simple  white  light  pioeessmg 
technique  hv  which  pseudocoloi  densitv  encoding  thiough  contrast 
reversal  can  be  easily  obtained  *"  He  showed  that  the  white-light 
pioeessmg  technique  otlcis  no  apparent  resolution  loss,  and  the 
system  is  veiv  simple,  versatile,  and  economical  to  opeiatc  Ihe 
while  light  pioeessmg  system  also  piov  ides  a  ditcct-v  lew  mg  capability 
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that  i)  essential  lor  practical  applications  lluwcvci.  the  technique 
still  sutlers  one  mu|or  drawback  it  is  not  a  ical-tiiiie  pseudocolor 
encoding  technique.  In  a  more  recent  article  Yu  described  a  lech 
niqueol  real-time  white-light  pseudocoloi  encoding  in  density  and  in 
spatial  (requeues  1,1 

We  will  now  describe  a  technique  ol  contrast  rcvci sal  minding  loi 
gras  level  photogiaphic  linages  We  assume  that  negative  and  positive 
image  transparencies  ol  the  same  object  are  available,  and  we  let  the 
encodings  take  place  on  a  Iresli  photogtapluc  him  bv  sequential 
recording  ol  the  negative  and  the  positive  image  transparencies  w  it  It  a 
Ronchi  grating  in  two  angular  positions  I  or  example,  the  lust  re¬ 
cording  is  exposed  with  the  negative  transparence  at  II  .  and  the 
second  recording  is  made  with  the  positive  transparence  at  a  90 

If  we  place  the  spat  la  lie  encoded  transparence  at  the  input  plane 
I’ |  of  a  white-light  processor  in  l  ig  I.  but  without  the  ditliuction 
grating,  the  smeared  I  ourier  spectra  are  dispersed  into  the  typical 
rainbow  colors  in  the  Courier  plane.  I\ 

In  pseudocolor  encoding,  we  Idler /no  lirst-oidei  smeaied  spec¬ 
tra  through  green  and  led  color  lillcis.  respectively.  and  the  image 
irradiance  at  the  output  plane  I’,  will  be 

l(x.y)  -  l,^(x.v)  t  Ip'(x.v)  .  ('ll 

where  ln^  and  Ip,  are  the  io/or  //c.gumv  and  re, I  ,  ,/lnr  /umiis  c 
image  irradiances  1  he  output  image  is  the  siipei position  ol  a  green 
negative  image  with  a  red  positive  image.  1  Inis  a  hioad  range  ot 
pseudocolor  density  images  can  be  obtained  hv  the  while-light  pio- 
cessing  technique.  I  he  selection  ol  color  filters  is  at  binary.  that  is. 
for  different  color  filters  one  would  obtain  dillcient  shades  ol 
pseudocolor-coded  images. 

In  experimental  demonstration  we  show  a  pseudocolor-encoded 
x-ray  picture,  as  in  I  Tg  14.  We  see  that  a  broad  range  ol  pseudocolors 
can  be  obtained  by  this  technique,  and  the  color  coded  image  appears 
tree  trom  coherent  artdact  noise  W  e  would  also  note  that,  in  pseudo- 
coloring.  multicolor  filters  can  be  implemented  in  the  (  ourier  spec¬ 
tral  bands,  l  or  example,  to  obtain  a  dillercnt  shade  of  pseudo¬ 
color.  one  may  insert  a  blue  filter  (in  blue  region)  w  it h  a  green  idler  (in 
green  region)  in  the  same  order  spectral  band  I  bus  a  In  oad  range  ol 
color  combinations  can  be  obtained  hv  this  pseudocoloi  mg 
technique. 

9.  conci.iidinc;  rim.yrks 

We  have  reviewed  the  basic  principles  ol  colicicnt  and  incoherent 
optical  processing  techniques.  We  have  shown  that  both  the  coherent 
and  incoherent  optical  processing  operations  are  able  to  be  analyzed 
with  linear  system  theory.  Although  coherent  optical  processing 
techniques  have  been  used  historically  for  most  ol  the  optical  infor¬ 
mation  processing  operations,  the  coherent  processing  system  is 
plagued  with  artdact  noise  which  frequently  limits  its  processing 
capability.  In  this  paper,  we  have  shown  th.it  llicic  aie  several  pio- 
cessuig  techniques  which  can  be  operated  with  incoherent  oi  white- 
light  source  Since  white-light  contains  all  the  v  isiblc  wavelengths,  it 
is  partieuliii ly  suitable  loi  coloi  photograph.  .mage  processing  In 
addition. the  white-light  source  is  generally  Iv.s  expensive  and  the 
system  stability  is  not  as  critical  as  with  a  coherent  processor. 

We  have  described  briefly  the  basic  concepts  ol  image  deblurring, 
image  subtraction,  nonlinear  processing,  space-vaiiant  pioccssing. 
color  image  retrieval  and  enhancement,  anil  pseudocolor  encoding 
ol  photogtapluc  images  with  column  oi  white  light  illiiuiuiulioii. 
We  stress  that  the  optical  processing  opeiation.  cither  coheicnt  oi 
incoherent,  is  primarily  based  on  the  I  our  icr  Irausloi  m  properties  ot 
lenses.  We  note  that  opt  it  al  plot  evunp  pc  ne  i  ally  opeiales  tu  a  lixetl 
parameter  mode,  and  lacks  the  llexibdily  ol  a  digital  conipulci 
Nevertheless,  there  still  are  processing  operations,  a  lew  ol  v\  Inch  we 
have  dlustialcd.  which  can  lie  achieved  with  analog  optical  leili 
iliques.  particularly  with  photographic  images  Since  spatial  resolu 
turn  is  very  unportaiil  in  photographic  iinagciy,  the  new  (lends  u. 
white-light  processing  can  be  expected  to  open  new  domains  in  color 
photographic  processing 


In  spite  ot  the  tlexibdity  ot  digital  image  processing,  optical 
methods  otlci  the  advantages  ol  grcatei  capacity,  simplicity,  and 
lower  oi.l  Instead  ol  commuting  each  other,  we  can  expect  a 
gi.ntu.il  niciguig  <4  the  optical  and  digital  techniques  1  he  continued 
dev  clop  me  in  ot  opr  ic.il -digital  inter  laces  and  input  dev  ices  will  lead 
lo  a  l  r  it  it  l  ul  i  suit  ol  livbrid  optical  digital  processing  techniques, 
iilili/uig  i he  siiengllis  ul  both  processing  operations 
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Source  encoding  for  image  subtraction 
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A  technique  nl  soum*  em-iMfing  <»(  an  extended  incoherent  source  f*»r  *uMr;u  lion  presented  The  sourer 

encoding  constraints  arc  obtained  troin  l  be  i  >  therein  v  requirement  lor  the  >ul»l  raet  i<  *n  *  »perai  ion  Sour*  e  encoding 
increases  the  available  light  power  tor  the  processing  operation,  as  a  small  nuoherent  source  is  no  longer  required 
Ail  experimental  result  obtained  wit  h  I  lie  encoded  incoherent  sour*  e  to*  tmique  is  given  A  similar  result  obtained 
by  using  a  coherent  technique  is  included  lor  comparison. 


Optical  image  subtraction  with  complex  amplitude 
was  described  by  (labor  <7  at  1  more  than  a  decade  ago. 
The  technique  involves  successive  recordings  of  two  or 
more  complex  diffraction  patterns  on  a  holographic 
plate  and  the  subsequent  reproduction  of  the  composite 
hologram  images.  A  few  years  later,  Bromley  <7  at  - 
described  a  holographic  Fourier  subtraction  technique 
with  which  a  real-time  image  and  a  previously  recorded 
hologram  image  can  he  subtracted.  Although  Bromley 
ft  at.  achieved  good  image  subtraction  in  their  experi¬ 
ments  it  appears  that  the  illumination  lor  the  hologram 
image  reconstruction  must  he  arranged  carefully.  In 
a  more  recent  paper,  Lee  ft  al.  {  proposed  a  method 
whereby  image  subtraction  and  addition  can  also  he 
achieved  by  a  diffraction-grating  technique.  This 
technique  involves  the  insertion  of  a  diffraction  grating 
in  the  spatial-frequency  domain  of  the  coherent  optical 
processor. 

However,  most  of  the  image  subtraction  techniques' 
require  a  coherent  source.  Such  sources  introduce  co¬ 
herent  artifact  noise  which  limits  their  processing 
capabilities.  In  previous  papersr’>,i  we  have  proposed 
a  technique  of  image  subtraction  that  requires  an  in¬ 
coherent  point  source.  However,  a  small  incoherent 
source  is  difficult  to  obtain  in  practice.  Nevertheless, 
this  difficulty  can  he  removed  with  the  source-encoding 
technique  that  is  discussed  in  this  Letter. 

Optical  image  subt  raet  ion  with  the  diffraction-grating 
technique  developed  by  Lee  <7  at  1  is  basically  a  one- 
dimensional  processing  operation.  Instead  of  using  a 
point  source  of  light,  one  can  use  a  line  source  perpen¬ 
dicular  to  I  lie  separa  I  ion  of  I  he  I  wo  input  object  trims 
pnreneies.  Since  the  image  subtraction  operates  upon 
the  corresponding  image  points  to  he  subl rai  led,  a 
strictly  broad  spatial-coherence  requirement  is  not 
needed.  Thus  it  is  possible  to  encode  an  extended 
source  in  order  to  obtain  a  reduced  point  -pair  spat  ial 
coherence  for  the  image  sulit ruction  operation. 

In  evaluating  the  spatial  coherence  requirement  for 
subt  ruction,  we  apply  partially  coherent  imaging  theory7 
ill  I  he  spill  lnl  IrequeiK  v  plane  /  '  ,  ol  nil  opl  leal  pioccaaoi , 
ns  shown  in  Fig.  I .  The  mill  uni  coherence  liind  ion  is 

Al;l(.V:i,  v'a>  =  .1.1  P.’fs..,  V '•.>)/(.*  ..)/* 

X  (.r'  dK/l.r  *;()/< *.,< x'.,,  t:i)il  v.'d.v'x,  ( I ) 


where  the  integration  is  over  the  inpu'  plane  a  j  and 
a  ,i  are  the  spat  ial -coordinate  syslei  is  of  /'  ,  and 
/iM.Vj.  a  I  is  the  complex  coherent  e  function  at  the 
input  plant  /'., /Ia.i  is  tile  input  oh  ect  functional 
and  can  lie  expressed  as/(.Vjl  =  Op  —  ht,i  +  It  .ia  ,  + 
hut;  ( *|(.v _■)  and  < J|t.v-d  are  the  two  input  object  trans¬ 
parencies.  and 


K jl  v ,v;|l  =  exp  1/2 ir  -'0  j  (2( 

is  the  transmittance  function  between  planes  IM  and  1  ‘ 

X  is  the  wavelength  of  the  light  source,  and  /  is  the  focal 
length  of  the  transform  lens  L_>. 

Thus  the  mutual  coherence  function  immediately 
after  the  diffraction  grating  (J,  with  a  spatial  period  d 
=  (X/l//i„,  can  he  shown  to  he 


Adl-V; i,  A  ';,)  = 


I  O  ^0 

exp  |(2tt  —x-a 


exp 


exp 


■'-'r  *  i|  - 


'-"A/'  ' 


;.)  h  II 


(HI 


where  we  ignored  the  dc  term  of  the  diffraction  grating 
and  (t  =  1/2  1 1  I  <  os|27r(/i,i/X/)A  |||  is  a  cosine  grating. 
We  note  that  for  complex  image  addition  of  sine  grating 
should  he  used.  The  image  intensity  at  the  output 
plane  l’i  is 


/(.V  |) 


|  f  fi  i(V;i,  .v'alexp  |/2ir  '  1  ‘a  ,|di  ,dv  „ 


( I) 


where  I  lie  iulegrat  ion  is  over  I  he  spatial  I'reqiieiicv  plane 
and  a  |  is  I  he  nut  pul  spat  ial  coordinate  system. 

Let  us  mbslUiUe  Kq.s.  (2)  and  l.'l)  into  Kq.  t  it  .and 
integrate  over  the  spat  ial -frequency  plane.  Considering 
only  the  image  terms  around  I  he  origin  of  the  output 
plane  I  *  (,  we  have 


/„(  a  ,1  -  |u  .('..,b„l|  |0,<  v  ,)  -  OAx  ,)| 

Ml  |,i  .r.-h„)|||c,(,  ,i| '  i  \<>.u  ,i|  '| 

(hi 


From  Kq  (M  we  see  that  the  first  lenn  is  pmpml  ional 
to  the  intensity  of  (lie  subtracted  image  and  the  second 
term  is  proport  ional  to  the  sum  ol  the  image  ir  radial  ices. 


Koprinled  from  Optic*  l<ettcr*,  Vol.  f>\  pnge  4;V2.  September,  l!»HI 
f  opvrigfd  I!IH|  liv  die  Opl ir.il  •rict y  n|  Anti- rit  h  hhI  ri'prini"'l  lt\  pi  r  •  if  iP<  enpvrii'lil  mwiih 


September  HIHI  /  Vnl.  (i.  No.  (1  /  OPTICS  LETTERS  453 


)<> 


where  |p(2/io)|  is  the  degree  of  spatial  coherence.  If  the 
decree  of  coherence  |p(‘2/ii)l|  is  high,  i  e.,  if  p(2h,,)  1, 

then  Eg.  (5)  reduces  to 


/o(-x4)  ~  |0,<x4)  -  (>.<x4)|-.  |/r,(2/i„)|  ^  1.  (til 


Thus  we  see  that  spatial  coherence  is  required  for 
every  pair  of  subtracted  image  points.  In  other  words, 
only  a  point-pair  spatial -coherence  requirement  j 
-  x'2)|  ^  | for  |x2  -  x'2|  =  2/q(  >s  needed  for  the  sub¬ 
traction  operation. 

In  source-encoding,  we  let  the  intensity  transmittance 
of  the  encoding  mask  he 


S(x,) 


,v 

£  rect 

n  *  I 


(7) 


a  multiple-slit  source,  where  jV  is  the  number  of  encoded 
slits,  s  is  the  slit  width,  and  <1  is  the  spacing  between  slits 
of  the  encoding  mask.  We  note  that  d  is  also  the  spatial 
period  of  the  grating  (1.  At  the  input  plane  1C,  the 
spatial  coherence  function  can  he  shown," 


p2(x2  -  x'2l 


sin 


x  >  -  x'2\ 
h„  ) 


N  sin 


X  sine 


f  T,  (X-2  -  X'i) 

dhn 


(8) 


where  d  =  \//h().  Prom  this  equation  we  see  that  the 
last  sine  factor  is  identical  with  the  single-slit  case, 
which  represents  a  broad  spread  of  coherence  over  (x2 
—  x'2).  However,  the  first  factor,  lor  large  values  of  N, 
converges  to  a  sequence  of  narrow  pulses.  The  locations 
of  the  pidses  (i.e.,  the  peaks)  occur  at  every  x  =  x2  —  x'2 
=  n  X//<h  which  yields  a  spatial-coherence  discrimina¬ 
tion  of  at \[/d)  over  the  input  plane  1C.  Thus  the 
multislit  source  encoding  not  only  provides  the  point- 
pair  coherence  needed  for  image  subtraction  but  also 
provides  a  higher  available  light  power  for  the  operation. 
In  other  words,  the  multislit  encoding  utilizes  the  light 
source  more  effectively  so  that  the  inherent  difficulty 
of  acquiring  a  small  incoherent  source  can  be  re¬ 
moved. 

In  our  experiment  a  mercury-arc  lamp  with  a  green 
filter  was  used  as  an  extended  incoherent  source.  A 
multislit  mask  was  used  to  encode  the  light  source.  The 
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Pig.  2.  Image  subtraction,  (a),  (ht  Input  object  transpar¬ 
encies,  (c)  subtracted  image  obtained  with  incoherent  tech¬ 
nique,  (d)  subtracted  image  obtained  with  coherent  tech¬ 
nique. 


slit  width  s  was  2.5  pm,  the  spacing  between  slits  was  25 
pm,  and  the  overall  size  of  the  mask,  which  contained 
about  100  slils,  was  about  2.5  mm  X  2.5  mm.  The  focal 
lengths  of  the  transform  lenses  were  200  mm.  A  liquid 
gate,  which  contained  the  two  object  transparencies  of 
size  6  mm  X  8  mm  was  inserted  immediately  behind  the 
collimator.  A  sinusoidal  phase  grating  with  a  period  of 
25pm  was  used  in  the  spatial-frequency  plane  P2.  The 
reparation  between  the  two  input  images  to  P2  was  12.2 
mm. 

In  our  experiments,  a  set  of  binary  images  as  shown 
in  Figs.  2(a)  and  2(h)  was  used  as  input  objects.  Figure 
2(c)  shows  the  subtracted  image  obtained  with  this 
source-encoding  technique.  For  comparison  the  sub¬ 
tracted  image  obtained  with  the  conventional  coherent 
processing  technique  is  shown  in  Fig.  2(d).  As  can  be 
seen,  we  obtained  better  artifact-noise  suppression  by 
using  the  incoherent  technique,  which  results  in  a  better 
subtracted  image. 

We  have  introduced  a  source-encoding  technique  for 
image  subtraction.  We  stress  that  the  concept  of  source 
encoding  may  he  extended  to  other  information  pro¬ 
cessing  operations.  We  also  note  that  image  subtrac¬ 
tion  with  the  encoded  extended  incoherent -source 
technique  is  generally  simple,  versatile,  and  economical 
to  operate.  It  may  offer  a  wide  range  of  practical  ap¬ 
plications.  In  addition,  the  technique  is  also  capable 
of  operating  in  a  real-time  mode. 
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Image  subtraction  with  an  encoded  extended 
incoherent  source 

S.  T.  Wu  and  F.  T.  S.  Yu 


A  technique  of  encoding  an  extended  incoherent  source  for  image  subtraction  is  presented.  The  source  en¬ 
coding  is  obtained  from  the  coherence  requirement  for  image  subtraction  operation.  Since  the  coherence 
requirement  is  a  point  pair  concept  for  image  subtraction  the  encoding  can  take  place  by  spatial  sampling 
an  extended  incoherent  source  with  narrow  slit  apertures.  The  basic  advantage  of  the  source  encoding  is 
to  increase  the  available  light  power  for  the  processing  operation,  so  that  the  inherent  difficulty  of  obtaining 
a  very  small  incoherent  source  can  be  alleviated.  Experimental  results  obtained  with  this  encoded  incoher¬ 
ent  source  are  given.  Comparisons  with  the  results  obtained  by  processing  technique  are  also  provided. 


I,  Introduction 

One  most  interesting  and  important  application  of 
optical  information  processing  must  he  image  sub¬ 
traction.  The  applications  may  he  of  value  in  urban 
development,  earth  resource  studies,  meteorology, 
highway  planning,  land  use,  inspection,  automatic 
tracking  and  surveillance,  etc.  Optical  image  sub¬ 
traction  may  also  apply  to  electrical  and  video  com¬ 
munications  as  a  means  of  bandwidth  compression. 
For  example,  it  is  only  necessary  to  transmit  the  dif¬ 
ferences  between  the  code  words  or  images  in  successive 
cycles  rather  than  the  whole  code  word  or  the  entire 
image  in  each  cycle  (e.g.,  TV). 

In  1965,  optical  image  synthesis  by  complex  ampli¬ 
tude  subtraction  was  first  described  by  Gabor  et  al. 1 
The  technique  involves  successive  recordings  of  two  or 
more  complex  diffraction  patterns  on  a  holographic 
plate  and  the  subsequent  reproduction  of  the  composite 
hologram  images.  A  few  years  later,  Bromley  et  a/.2 
described  a  holographic  Fourier  subtraction  technique, 
for  which  a  real-time  image  and  a  previously  recorded 
hologram  image  can  he  subtracted.  Although  good 
image  subtraction  by  their  experiments  had  been  re¬ 
ported,  it  appears  that  the  illumination  for  the  hologram 
image  reconstruction  must  be  arranged  carefully.  In 
1970,  Lee  et  al.  '  proposed  a  technique  so  that  image 
subtraction  and  addition  can  also  he  achieved  by  a  dif¬ 
fraction  grating  technique.  This  technique  involves 
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insertion  of  a  diffraction  grating  in  the  spatial  frequency 
domain  of  the  coherent  optical  processor.  Good  results 
by  the  diffraction  grating  technique  were  also  reported 
in  their  article.  In  a  more  recent  article.  Zhao  et  at  * 
also  proposed  a  technique  of  image  subtraction  utilizing 
a  halftone  screen  method.  However,  thJr  technique 
cannot  he  implemented  in  real  time. 

There  are  several  other  techniques  available  for  image 
subtraction  which  can  he  found  in  a  review  paper  by 
Ebersole.5  However,  most  of  the  optical  information 
processing  techniques  require  a  coherent  source  to  carry 
out  the  subtraction  operation.  But  coherent  optical 
processing  systems  are  plagued  with  coherent  artifact 
noise,  which  frequently  limits  their  processing  capa¬ 
bilities. 

We  have  in  previous  papers*’-7  proposed  a  technique 
of  optical  processing  with  an  incoherent  source  for 
complex  signal  detection  and  image  deblurring.8-9  We 
have  also  extended  the  technique  for  possible  applica¬ 
tion  to  image  subtraction. ,on  However,  to  obtain  a 
spatial  coherence  requirement  for  subtraction  operation 
a  very  small  source  size  is  needed,  blit  a  small  incoherent 
source  is  difficult  to  obtain  in  practice.  Nevertheless, 
this  difficulty  may  he  alleviated  with  a  source  encoding 
(i.e.,  spatial  sampling)  technique,  so  that  the  extended 
incoherent  source  can  be  used.  We  have  briefly  dis¬ 
cussed  in  a  recent  communication  that  image  subtrac¬ 
tion  can  indeed  obtain  with  an  encoded  extended 
source.12  The  basic  object  ive  of  source  encoding  is  to 
utilize  the  light  power  more  effectively  so  that  image 
subtraction  can  he  carried  out  with  an  extended  inco¬ 
herent  source.  Moreover,  with  t  he  use  of  an  incoherent 
source,  coherent  artifact  noise  can  he  avoided.  We 
stress  that  Ibis  image  subtraction  system  is  capable  of 
operating  in  the  real-time  mode,  and  it  is  generally 
simple,  versatile,  and  economieal. 
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Kig.  1.  Young's  experiments  with  an  extended  source. 


Basically,  the  optical  image  subtraction  is  a  1-D 
processing  operation.  Instead  of  utilizing  a  point  source 
of  light,  a  line  source  of  light  can  he  used  for  the  sub¬ 
traction  operation.  Since  the  spatial  coherence  re¬ 
quirement  for  subtraction  operation  is  a  point-pair 
concept,  a  strictly  coherence  requirement  is  not  needed. 
In  other  words,  it  is  possible  to  encode  an  extended  in¬ 
coherent  source  to  obtain  a  point  -pair  coherence  re¬ 
quirement  for  image  subtraction  operation. 

II.  Young's  Experiment  with  Extended  Source 

We  will  now  illustrate  the  concept  of  Young’s  exper¬ 
iment  for  the  source  encoding.  Let  us  consider  first  a 
narrow  slit  of  light  source  .S,  situated  in  plane  P\  to  Fig. 

1 .  To  maintain  a  high  degree  of  spatial  coherence  be¬ 
tween  apertures  Qi  and  Qj,  the  source  £>i  should  be  very 
narrow.  In  other  words,  if  the  separation  between  the 
two  open  apertures  is  larger,  the  narrower  incoherent 
source  S i  is  required.  It  can  he  showm  that,  to  maintain 
a  high  degree  of  coherence,  the  slit  size  can  be  approx¬ 
imated  by1''*'14 

S  *  1X/12M  R,  (l) 

where  R  is  the  distance  between  planes  P\  and  P 2. 

We  now  consider  two  narrow  slits  of  light  sources  S i 
and  S*2,  as  shown  in  Fig.  1.  If  the  separation  d’between 
the  two  sources  Si  and  S>  satisfies  the  relation 

r\  ~  r'l  a  (ri  -  r-j)  +  m A,  (2) 

where  the  r’s  are  the  distances  from  sources  .S',  and  S». 
to  the  open  apertures  and  us  shown  in  the  figure, 
m  is  an  arbitrary  integer,  and  \  is  the  wavelength  of  the 
light  source,  the  interference  fringes  due  to  each  slit 
source  are  in-phased,  and  u  brighter  fringe  pat  tern  can 
be  observed  at  plane  /’a.  With  the  application  of  Kq. 
(2),  we  can  employ  many  narrow  slit  sources  since  we 
wish  to  obtain  a  coherent  fringe  pattern  at  the  output 
plane  P,\.  We  note  that  the  separations  between  any 
of  the  two  slit  sources  should  satisfy  tin  fringe  (i.e., 
spatial  coherence)  condition  of  Kq.  (2).  If  the  separa¬ 
tion  K  lietween  planes  / *,  and  is  large,  i.e.,  R  W  </,  /aid 
R  »  2/i„,  the  coherent  condition  of  Kq.  (2)  becomes 

./•  =  «r.  |,\H/l2»i„»l.  C») 


From  this  equal  ion  we  see  that  equal  spacing  slits  can 
he  used  so  that  a  brighter  fringe  pattern  can  be  ob¬ 
served.  We  note  that  the  intensity  of  the  fringe  pattern 
increases  linearly  as  the  number  of  slits  increases.  Thus 
on  one  hand  the  source  encoding  preserves  the  coher¬ 
ence  requirement,  and  on  the  other  hand  it  increases  the 
overall  intensity  of  illumination.  Therefore,  with  ap¬ 
propriate  source  encoding,  an  extended  source  may  be 
efficiently  utilized. 


III.  Spatial  Coherence  Requirement 

W  e  will  now  adopt  the  concept  of  source  encoding  in 
evaluating  the  spatial  coherence  requirement  for  image 
subtraction  operation.  With  reference  to  the  incoher¬ 
ent  optical  processor  of  Fig.  2,  we  see  that  the  processor 
is  similar  to  that  of  a  coherent  optical  processor  except 
with  an  extended  incoherent  source  and  an  encoding 
mask.  Since  image  subtraction  is  a  1  -D  operation,  we 
will  adopt  a  1-1)  notation  for  our  analysis. 

In  evaluating  the  spatial  coherence  requirement,  we 
use  partially  coherent  imaging  theory.1-'14  The  mutual 
coherence  function  at  the  spatial  frequency  plane  P.\ 
is 

teilx.i.x  )  =  j*J  M  J  a/U-d/'U'j)K  jU2,x,|t 

x  K[,  u,,.v  itd.tjdx'  .,  (41 

where  the  integration  is  over  the  input  plane  P>\  .v-j,  x->. 

and  .t;,  are  the  position  coordinates  of  P±  and 
respectively,  gjl.vj.v)  is  the  complex  coherence  function 
at  the  input  plane  /*•_>;  f t-v •_»>  is  the  input  function  at 
which  can  be  expressed  as 

fix-i )  =  (tiix-y  ~  /| 0>  +  0  >(X‘l  +  (M 

where  0\(x>)  and  Oi(-Vo)  are  the  two  input  object 
transparencies  and 


K  al  x  2-x  ,\) 


U>l 


is  the  transmittance  function  between  planes  l'<  and  P.\, 
\  is  the  wavelength  of  the  light  source,  and  /  is  the  focal 
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length  of  the  transform  lens  Equation  (-1)  can  also 
he  written 


,1  =  J  |  p..(A  ■_>,>  _l|(  -  ll,,)  +  <>J  (X_.  +  (lot] 

x  (nj  u.j  -  /tll)  +  it]  U_,  +  /i„i| 

(A  l-V  t  -  a  ,.x  \ 

I  ‘J  IT  — — - jit*  •_.</.»  til 

where  superscript  *  denotes  the  complex  conjugate. 

It  is  clear  that  the  mutual  coherence  function  im¬ 
mediately  behind  the  diffraction  granting  (1,  with  a 
spacing  period  d  =  (\/)//i,i,  is 

M  ,(A  i.a  ,1  «  jexp|i  2ir  ~  a  , |  -  exp|-i'2r  ^  x;tJ| 

X  exp|— (gx  x  ,j  -  cxp|igir  *  >j  M: i(x.i.Xj),  (HI 

where  we  ignore  the  dc  term  of  the  diffraction  grating. 
The  image  intensity  at  the  output  plane  /*4  is 

/up  =  Jj  (J  ,(X:i,Xj>  »'X|]|rJrr  — — — -  x  ,Ji/ai|/a  ,,  (91 

where  the  integration  is  over  the  spatial  frequency 
plane,  and  ,v4  is  the  output  spatial  coordinate  system. 

By  substituting  K.qs.  (7)  and  (8>  into  Kq.  t!))  and  in¬ 
tegrating  over  the  spatial  frequency  plane,  we  have 

llxtl  =  ptx  >,A2)|0||A3  -  ha)  +  Oj(aj  +  /iiil||t/|Us  —  ho)  +  0,x,  - 

•  | it ( A 2  E  X4  +  /l„)fltX2  +  A 4  +  ho)  +  JlXj  +  A4  —  /|0)6(Ai  +  A 4  -  / 
-  it  (A  -j  +  A  4  +  /l„)lt(A,  +  A4  -  llo)  it  ( A  -2  +  A4  -  /|||)6(A2  +  A  4  +  / 

where  (5(x)  is  the  Dirac  delta  function.  Let  us  assume 
that  n(x>,x>)  takes  the  form  p(.t-»  -  x>)  and  gtil  = 

If  we  evaluate  Kq.  ( 10)  termwise  and  note  that 
the  input  images  are  spatially  limited,  we  have 

/(A4»  =  ptl/llbM  — X4)|-  +  |fM-xO|'-!  -  p ( *i/l ut( 1 1 1  X  4  )02  (-A4) 

-ft’  ( ’2h)))( )|(  — A 4)/ )2t  A 4)  +  p(l>)|()|(-X4  -  2/l | - 

+  1<M-a,  +  2M-I,  tin 

where  p(2/to)  =  |ji(2/io)|  exp (/</>)  is  a  complex  quantity. 
We  stress  that  phase  factor  </>  can  he  avoided  by  ad¬ 
justing  the  grating  position  (». 

We  now  consider  only  the  image  terms  around  the 
origin  of  the  output  plane  P 4, 

/|>(  A 4 )  =  |Mt2*„>||(),(A«)  -0...(A4ir- 

+  |1  -  |mI2/i11»1|||0|(a.,»|-  +  |(>2(x.|)|  j.  for,/.  =  0. 


From  Kq.  (12)  we  see  that  the  first  term  is  propor¬ 
tional  to  the  intensity  of  the  subtracted  image,  and  the 
second  term  is  proportional  to  the  sum  of  the  image  ir- 
radiances,  where  |p(2hp)|  is  the  degree  of  spatial  co¬ 
herence.  If  the  degree  of  coherence  |p(2/i„)|  ishigh.i.e., 
n(2h{t)  1,  Kq.  (12)  reduces  to 

/,it-A4*~  |f  >i(a.,  )-<>■.(•»,  M-.  for  |/iC-’/i„)|  «  t.  tt:n 

'I'hus  wi>  see  that  the  spatial  coherence  is  only  needed 
for  every  pair  of  points  |or ■,  -  x,|  -  2/»t).  In  other  words, 
only  a  point  pair  spatial  coherence  is  required  for 
suht  rad  ion  operation. 


IV.  Source  Encoding 

We  will  now  search  a  source  encoding  so  that  point- 
pair  spatial  coherence  can  he  established.  We  will 
adopt  the  concept  of  Young’s  experiment  that  we  de¬ 
scribed  in  a  previous  section.  Now  we  insert  a  mask 
transparency  lor  the  source  encoding  at  the  front  local 
plane  l\  of  the  collimator  /.  1  as  show  n  in  Fig.  2.  The 
spatial  coherent  function  p ( v . I  over  the  input  plane 
can  he  written1 1 

fill  ...A  .  I  =  I  .S'lA  iIK"|IA  ,.a  2iK|(X|.X2u/X|.  I  141 

where  Nl.t  i)  is  the  intensity  transmittance  function  ot 
the  mask,  and  K ,(.v  ,,x-j)  is  the  transmittance  function 
between  planes  P\  and  P<.  We  assume  that  the  mask 
is  located  within  an  isoplanatic  patch,  and  K  s (.v  j ,.v_>)  can 
be  written111 


K  |(A  |.A_i)  =  t-xpi  2ir  ■  +  t  j  A —  A  |  — )]  •  (IS) 

where  t(x)  is  the  wave  aberration  of  the  collimator,  and 
b  is  the  distance  between  the  collimator  and  the  input 
plane  P  >■ 

We  note  that  if  6  is  sufficiently  small,  i.e.,  (b/f)  *  Imax 
«  a 2,  the  transmittance  function  of  Eq.  (15)  can  reduce 
to 


h0) 

/|,iI|(/a'2</.V  -j.  I  1(1) 


K  |(a  i.aj)  expji  |'2tt  —j:  +  ( (  a  2 )  j  I  ■  lid) 

By  substituting  Eq.  (16)  into  Eq.  (14),  the  spatial  co¬ 
herence  function  becomes 

p(A  •„•,.(  2)  *S  4*Xpi  1< (A'2l  -  <lx  >)| 

X  J  .S'(ai)  exp  (2ir  ^  (xs  -  x2)Jdxt  (17) 

From  the  above  equation  we  see  that  the  spatial  co¬ 
herence  function  is  the  Fourier  transform  of  the  mask 
transmittance  function  modulated  by  a  phase  (wave 
aberration)  factor.  However,  the  phase  aberration  will 
not  affect  the  degree  of  mutual  coherence  |^i(.Vo,.va)  l-  In 
the  remaining  analysis,  we  shall  ignore  the  phase  aber¬ 
ration  and  assume  that  ix{x>,x2)  takes  the  form  p(.V2  ~ 
*•;). 

Now  we  evaluate  the  degree  of  coherence  |p(2/to)|  for 
two  different  cases.  We  will  first  evaluate  a  single-slit 
encoding,  i.e., 

S(.t|)  =  rect  (a i/s),  (18) 

where  s  is  the  slit  width.  By  substituting  Kq.  (18)  into 
Kq.  (17),  we  have 

plA-.i  -  A  j )  =  Sim  — (A  .1  -  A  .i),  (19) 

where  the  phase  factor  was  ignored.  Since  the  spacing 
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Table  I.  Spatial  Coherence  Requirement  lor  Slngle-Slil  Math 

,  'il  |  /'J  I  /:'»  I  III  I  "-'ii 

pfjl  t„|  (I  (I.V.Mi  n  tetri  ilium 


period  of  the  unit in^  (».d  =  t /  \  > //i , >.  if  we  let  x2  —  x2  = 
'2h, i  =  (2/A)  '</.  Ki|.  I 19)  t  an  lit'  written 


Thus  we  see  that  the  decree  of  spatial  eoherenee 
| 1  -/l 1 1 1 1  ttepeiuls  upon  the  ratio  of  the  slit  width  s  to  tlie 
spacing  d. 

rI\i  gain  a  feeling  of  magnitude,  we  provide  several 
values  of  pt2/t,>)  in  Table  I,  from  which  we  c  n  see  that 
a  high  degree  of  spatial  coherence  can  he  attained  only 
through  a  very  narrow  slit.  For  example,  it  the  spacing 
of  the  grating  d  =  25  pm,  to  achieve  a  high  degree  of 
spatial  coherence  a  slit  width  s  <  2.5  /.on  should  lie  used. 
Thus  it  makes  the  source  too  weak  lor  a  practical  pro¬ 
cessing  operation. 

As  noted  in  the  previous  section,  the  spatial  coherence 
requirement  for  image  subtraction  is  a  point  pair 
problem.  It  is  possible  to  encode  the  extended  source 
with  number  of  narrow  slits.  Thus  with  a  mult islit 
source  encoding,  an  .V  fold  light  power  can  he  used  for 
the  image  subtraction  operation. 

We  will  now  let  the  intensity  transmittance  of  the 
encoding  mask  he 

;v  l\  -  n<r\ 

.Step  =  recti — - - 1  •  l-'tl 

where  s  is  the  slit  width,  and  ti‘  is  the  spacing  between 
slits. 

Hy  substituting  Kq.  (21)  into  Kq.  (17),  the  spatial 
coherence  function  becomes 


where  x  =  x 2  —  From  the  above  equation  we  see 
that  the  last  sine  factor  is  identical  to  the  single-slit  case 


of  Kq.  (19),  which  represents  a  broad  spread  of  coher¬ 
ence  over  x.  However,  the  first  factor  for  large  values 
of  N  converges  to  a  sequence  of  narrow  pulses.  The 


locations  of  the  pulses  (i.e.,  the  peaks)  occur  at  every  .r 
=  *2  ~  x-2  *  n(\f/d‘).  Thus  this  factor  yields  the  fine 
spatial  coherence  discrimination  at  every  point-pair 
separated  at  distance  (\f/<i’)  over  the  input  plane  I’>. 

If  we  let  the  spacing  of  d'  equal  the  spacing  d  of  the 
diffraction  grating  (1  (i.e.,  d'  =  d),  the  spatial  coherence 


of  Kq.  (22)  becomes 


where  we  substitute  d  =  ( X/)//i«.  From  Kq.  (2:t),  we  see 
that  a  sequence  of  narrow  pulses  occurs  at  x  =  x-<  —  x2 
=  nh 0,  where  11  is  tin  integer,  and  their  peak  values  are 


Kin  :t.  Coherence  ftii«' lion  uhiainetl  with  mult  islit  source  encoding, 
whore  i  =  |  r_,  -  t  |  and  s/d  -  1.5 


weighed  hy  a  broader  sine  factor,  as  shown  in  Fig.  2.  It 
can  lie  shown  that  the  width  of  the  pulses  is  inversely 
proportional  to  the  number  of  slits  .V.  Thus  the  mul- 
I islit  source  encoding  not  only  provides  a  point  pair 
coherence  requirement  for  image  subtraction  but  also 
a  higher  available  light  power  for  the  operation.  In 
other  words,  the  multislit  encoding  utilizes  the  light 
source  more  efficiently,  so  that  the  inherent  difficulty 
of  acquiring  a  small  incoherent  source  can  be  alle¬ 
viated. 

V.  Temporal  Coherence  Requirement 

So  far  we  have  considered  only  the  quasi-mono- 
chromatic  light,  hut  the  effect  of  the  temporal  coherence 
has  not  been  discussed.  Since  the  scale  of  Fourier 
spectrum  varies  with  wavelength,  there  is  a  temporal 
coherence  requirement  for  every  processing  operation. 
With  this  consideration,  we  must  limit  the  temporal 
bandwidth  A  A  of  the  source  so  that  the  dispersed  Fou¬ 
rier  spectra  will  not  spread  beyond  the  allowable  limit. 
In  the  image  subtract  ion  operation,  we  should  limit  the 
spectrum  spread  within  a  very  small  fraction  of  the 
grating  spacing  d.  i.e., 

Iqin./AXtcjirl  «  d,  U24) 

where  p„,  is  the  highest  angular  spatial  frequency  of  the 
input  objects,  f  is  the  focal  length  of  the  transform  lens, 
and  AA  is  the  spatial  bandwidth  of  the  source. 
Therefore,  the  temporal  bandwidth  of  the  source  should 
he  limited  hy  the  following  inequality: 


where  A  is  I  be  center  wavelengt  h  of  the  light  source,  and 
2/ih  is  t  he  separat  ion  oft  lie  input  images. 

To  gain  a  practical  feeling,  we  let  /in  =  (>.(>  mm,  X  = 
r»-H»l  A  and  take  a  factor  of  10  of  Kq.  (2f>).  The  temporal 
bandwidth  requirements  AA  for  various  values  of  spatial 
frequencies  tire  tabulated  in  'Fable  11.  We  see  that, 
if  the  spatial  frequency  of  the  input  objects  is  low,  a 
broader  temporal  bandwidth  of  the  light  source  can  he 

I  l  M'  oihIh.i  lam  /  v,,l  ;■!>.  No  ;n  /  AH'I  II  I)  Ol’lll;:, 


4l)!l.'j 


ritf-  -I.  I  ritftf'r  jmhfrmlion;  mill  h  hum  w  Iinm*  objrrl;  b»>  mid  (l>)  injMil  olifcrl  trm»*|>mviu  ii>;  ii  \  »uhlnirltil  imu^r  nhbimril  wilh  mi  imolii-rritl 

h'l  llltlqlli',  ( ( 1 1  till  hi  |  m  |i<i|  HIM)  •  ithl  ill  lir<|  \\  1 1 1 1  .1  I  I  >||I  l  nil  In  ltnii|l|r, 
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1C) 

Fiji.  5.  Imosie  subtraction  parkin)!  lot:  tal.  (li)  inpm  object  transparencies;  <t  I  subtracted  image  obtained  with  an  incoherent  technique. 


this  figure,  the  profile  of  the  missing  passenger  car  can 
readily  he  seen  at  the  output  image  plane  It  is  also 
interesting  to  note  that  the  parking  line  on  the  right  side 
of  the  missing  passenger  ear  and  the  shadow  can  clearly 
he  seen  in  the  subtracted  image  of  Fig.  5(c). 

VII.  Summary 

We  have  introduced  a  source  encoding  technique  for 
image  subtraction.  The  source  encoded  function  was 
evaluated  from  a  specific  coherence  requirement  for 
image  subtraction  operation.  Since  the  image  sub¬ 
tract!.  ,  is  a  point-pair  problem,  and  the  processing  is 
essentially  a  1  -I)  operation,  it  is  possible  to  obtained  a 
required  coherence  function  by  encoding  an  extended 
source  with  a  set  of  narrow  slits.  In  other  words,  the  slit 
width,  the  spacing  of  the  slits,  and  the  number  of  slits 
determine  the  spatial  coherence  requirement. 

The  basic  advantage  of  the  source  encoding  is  to  in¬ 
crease  the  available  light  power  for  the  processing  op¬ 
eration,  so  that  the  inherent  difficulty  of  acquiring  a 
very  small  incoherent  source  can  be  alleviated.  We 
stress  that  the  concept  of  source  encoding  may  be  ex¬ 
tended  toother  incoherent  information  processing  op 
(Till  ions. 

Aside  from  the  spatial  coherence  requirement,  there 
is,  however,  a  lcnipcr.il  coherence  requirement  Ibr  the 


image  subtraction  operation.  If  the  spatial  frequency 
requirement  for  the  image  subtraction  is  high,  a  higher 
tempera!  coherence  (i.e.,  a  narrower  spectral  width)  of 
the  light  source  is  required. 

In  experimental  demonstrations,  we  have  shown  both 
the  results  obtained  with  the  incoherent  processing  and 
coherent  processing  techniques.  By  comparing  these 
results,  we  conclude  that  the  incoherent  processing 
technique  offers  a  better  artifact  noise  suppression  and 
better  image  quality. 

Finally,  we  also  stress  that  the  incoherent  image 
subtraction  technique  is  generally  very  simple,  versatile, 
and  economical  to  operate.  It  may  offer  a  wide  range 
of  practical  applications.  In  addition,  the  technique  is 
also  capable  of  operating  in  a  real-time  mode. 
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Abstract.  A  relation  between  spatial  coherence  function  and  source  encoding  intensity 
transmittance  function  is  presented.  Since  the  spatial  coherence  is  depending  upon  the 
information  processing  operation,  a  strictly  broad  spatial  coherence  function  may  not  be 
required  for  the  processing.  1  he  advantage  of  the  source  encoding  is  to  relax  the  constraints 
of  strict  coherence  requirement,  so  that  the  processing  operation  can  be  cai  i  led  out  with  an 
extended  incoherent  source.  I  mphasts  of  the  source  encodings  and  experimental  de¬ 
monstrations  are  given.  1  he  constraint  of  tempcral  coherence  requirement  is  also  discussed 

I’ACS:  42.30,  42.80 


I  he  use  of  coherent  light  enables  optical  processing 
systems  to  carry  out  many  sophisticated  information 
processing  operations  |  I.  2|  However,  coherent  opti¬ 
cal  processing  systems  are  contaminated  with  coherent 
artifact  noise,  which  frequently  limits  their  processing 
capabilities.  Recently,  attempts  of  using  an  incoherent 
source  to  carry  out  complex  information  processing 
operations  had  been  pursued  by  several  investigators 
[3  6],  The  basic  limitations  of  using  incoherent  source 
for  partially  coherent  processing  is  Ihe  extended  source 
si/e.  I  o  achieve  a  broad  spatial  coherence  function  at 
the  input  plane  of  an  optical  information  processor,  a 
very  small  source  si/e  is  required.  However,  such  a 
small  light  source  is  difficult  to  obtain  in  practice.  We 
have,  nevertheless,  shown  in  recent  published  papers 
1 7  |0]  that  there  are  information  processing  oper¬ 
ations  which  can  be  carried  out  with  incoherent 
source.  In  other  words,  a  strictly  broad  coherence 
requirement  may  not  be  needed  for  some  optical 
information  processing  operations. 

In  this  paper,  we  shall  describe  a  linear  transformation 
relationship  between  spatial  coherence  function  and 

*  ViMtiitg  scholar  fiom  Shanghai  (Iphcal  liisiiuincul  Kcscaich 
liistiiiitc.  Sluiigh.ii.  China 

"  ViMlmg  siholai  1 1 •  tin  Shanghai  Inslihtic  nl  0|»ui%  .nut  In. 
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source  encoding  intensity  transmittance  func'  on. 
Since  the  spatial  coherence  rcquiicmcnl  is  depending 
upon  the  information  processing  operation,  a  more 
relaxed  cohcience  function  may  he  used  for  a  specific 
processing  operation.  Uy  I  miner  transforming  this 
cohcience  lunctiou.  a  source  encoding  intensity  tians- 
miltance  function  may  be  found, 
ihe  purpose  of  source  encoding  is  to  reduce  the 
coherent  requirement,  so  that  an  extended  incoherent 
source  can  be  used  for  the  processing.  In  other  words, 
the  source  encoding  technique  is  capable  of  generating 
an  appropriate  coherence  function  for  a  specific  infor¬ 
mation  processing  operation  and  at  the  same  time  it 
utilizes  the  available  light  power  more  effectively.  We 
shall  illustrate  examples  that  complex  information 
processing  operation  can  actually  be  curried  out  by  an 
encoded  extended  incoherent  source.  I  xperimental 
illustrations  with  this  source  encoding  technique  are 
also  included. 

Source  I  ncmling  with  Spatial  Coherence 

We  shall  begin  um  discussion  with  the  Young's  experi¬ 
ment  under  extended  incoherent  source  illumination, 
as  shown  in  I  ig  I  l-'irst,  we  assume  that  a  narrow  slit 
is  pi, ned  at  plane  I’,  behind  an  extended  souice  lo 
maintain  a  high  degree  of  spatial  coherence  between 
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large,  i.e..  R  ./  aiul  /<>>2/i„.  llien  llie  spacing  </ 
between  1  he  souice  slits  becomes. 
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SOURCE  4  |S. 


fig  1  Youngs  experiment  with  exteiiilal  source  illtuniiuiiion 


the  slits  Q,  and  („),  at  I’,,  it  is  known  that  the  source 
si/e  should  be  very  narrow.  IT  the  separation  between 
Q,  and  (,),  is  large,  then  a  narrower  slit  si/e  S,  is 
required,  l'hus,  to  maintain  a  high  degree  til  spatial 
coherence  between  Q,  and  Q,.  the  slit  width  should  be 

[II] 


where  R  is  the  distance  between  planes  I’,  and  I’,,  and 
2h„  is  the  separation  between  (.),  and  tlig.  I). 
l  et  us  now  consider  two  narrow  slits  of  S,  and  S, 
located  in  source  plane  I*,.  We  assume  that  the 
separation  between  S,  and  S,  satisfied  the  following 
path  length  relation : 

ri  -  r\  =  <r,  -  r,l  + m/.,  (2) 

where  the  r's  are  the  respective  distances  front  S,  and 
S,  to  Q,  and  Q,.  as  shown  in  the  figure,  m  is  an 
arbitrary  integer,  and  /.  is  the  wavelength  of  the 
extended  source.  Then  the  interference  fringes  due  to 
each  of  the  two  source  slits  S,  and  S,  would  be  in 
phase.  A  brighter  fringe  pattern  can  be  seen  a!  plane 
l’3.  To  further  increase  the  intensity  of  the  fringe 
pattern,  one  would  simply  increase  the  number  of 
source  slits  in  appropriate  locations  in  the  source  plane 
P,  such  that  every  separation  between  slits  satisfied  the 
coherence  or  fringe  condition  of  (2).  If  separation  R  is 


Front  the  above  illustration,  we  see  that  by  propeily 
encoding  an  extended  source,  it  is  possible  to  maintain 
the  spatial  coherence  between  (./,  and  (J2.  and  at  the 
same  time  it  increases  the  intensity  of  illumination. 
Thus,  with  a  specific  source  encoding  technique  for  a 
given  information  processing  operation  may  result  a 
better  uttli/atioii  of  an  extended  source. 

To  encode  an  extended  source,  we  would  first  search 
for  a  spatial  coherence  function  for  an  information 
processing  operation.  With  reference  to  an  extended 
source  optical  processor  of  I  ig.  2,  the  spatial  coherence 
function  at  input  plane  P,  can  be  written  (11] 

/  lx,.  x,l  -  |  |  .S'fx,  \K  ,(\,.  \j)K,(x,.  N.li/x, .  (4) 

where  the  integiation  is  over  the  source  plane  P, .  .Six,  I 
is  the  intensity  transmittance  function  of  a  source 
encoding  mask,  and  /v,(\,.x2)  is  the  transmittance 
function  between  source  Plane  P,  the  input  plane  P,. 
which  can  be  written 

fC , lx , . x , 1 expj i ( 2ir X ! i  .  (5) 

By  substituting  /s,U|.x,|  into  (41.  we  have 

/'(x,-x’,)=  j  j  St X , ) ex p I i 2 ^ ( x J - x 2 1  </x, .  (6) 

From  the  above  equation,  we  see  that  the  spatial 
coherence  function  and  source  encoding  intensity 
transmittance  function  forms  a  Fourier  transform  pair 

si x j )  — . ^ \ ,  -  \ , t] .  (7) 

and 

/•(x,-x',)  =  >  ‘(six,)].  |8) 

where  .*  denotes  the  Fourier  transformation  oper¬ 
ation,  If  a  spatial  coherence  function  for  an  infor¬ 
mation  processing  operation  is  provided,  then  the 
source  encoding  intensity  transmittance  function  can 
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be  determined  through  Fourier  transformation  of  (7). 
We  note  that  the  source  encoding  function  N(x,(  can 
consist  of  apertures  or  slits  of  any  shape.  We  further 
note  that  in  practice  ,S'(x,)  should  be  a  positive  real 
function  which  satisfies  the  following  physical  realiz¬ 
able  condition : 

OgS(.\,)Sl.  (»■)) 

For  example,  if  a  spatial  coherence  function  for  an 
information  processing  operation  is 

n.v, -  Vj )  =  rect |'V  1  -  V^|,  (1 0) 

where  -I  is  an  arbitrary  positive  constant,  and 

l.vl  II.  \x\^A, 
ru  |  It),  otherwise. 

then  the  source  encoding  intensity  transmittance 
would  be 


Since  .V(v,|  is  a  bipolar  function,  therefore  it  is  not 
physically  realizable. 


Temporal  Coherence  Requirement 

There  is,  however,  a  temporal  coherence  requirement 
for  incoherent  source.  In  optical  information  process¬ 
ing  operation,  the  settle  of  the  Fourier  spectrum  varies 
with  wavelength  of  the  light  source,  therefore,  a 
temporal  coherence  requirement  should  be  imposed  on 
every  processing  operation.  If  we  restrict  the  Fourier 
spectra,  due  to  wavelength  spread,  within  u  small 
fraction  of  the  fringe  spacing  </  of  a  complex  spatial 
filter  (e.g„  deblurring  filter),  then  we  have. 


2  n 


(12) 


where  I  <f  is  the  highest  spatial  frequency  of  the  filler, 
P,„  is  the  angular  spatial  frequency  limit  of  the  input 
object  transparency,  f  is  the  focal  length  of  the  trans¬ 
form  lens,  and  ,17.  is  the  spectral  bandwidth  of  the  light 
source.  The  spectral  width  or  the  temporal  coherence 
requirement  of  the  light  source  is,  therefore. 


where  7  is  the  center  wavelength  of  the  (iplu  souice,  I’ll,, 
is  the  size  ol  the  input  oli|ccl  liauspaiency,  and 

2/t„  -  (7.  /  )/,/. 

In  order  to  gain  some  feeling  of  magnitude,  we  provide 
a  numerical  example,  l  et  us  assume  that  the  si/e  ol  the 
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object  is  2/i()  =  5  min.  the  w  avelength  of  the  light  source 
is  7.-54MA.  and  we  take  a  factor  10  for  (Hi  for 
consideration,  that  is 


.17.= 


I  On/. 
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Several 
various 
Table  I. 


values  of  spectral  width  requirement  .17.  for 
spatial  frequency  Vm  are  tabulated  in 


Front  f  able  I.  we  see  that,  if  the  spatial  frequency  of 
the  input  object  transparency  is  low.  a  broader  spectral 
width  of  light  source  can  he  used.  In  other  words,  if 
higher  spatial  frequency  is  required  for  an  information 
processing  operation,  then  a  narrower  spectral  width 
of  light  source  is  needed. 


(examples  of  Source  Fncoding 

We  shall  now  illustrate  examples  of  source  encoding 
for  partially  coherent  processing  operations.  We  would 
first  consider  the  correlation  detection  operation 

ri2i. 

In  correlation  detection,  the  spatial  coherence  require¬ 
ment  is  determined  by  the  size  of  the  detecting  object 
(i.e.,  signal).  To  insure  a  physically  realizable  encoded 
source  transmittance  function,  we  assume  a  spatial 
coherence  function  over  the  input  plane  I*,  is 


where  .1,  is  a  first-order  Bcsscll  function  of  first  kind, 
and  li„  is  the  size  of  the  detecting  signal.  A  sketch  of  the 
spatial  coherence  as  a  function  of  |  v ,  -  v',|  is  shown  in 
Fig.  3a.  By  taking  the  Fourier  transform  of  (15).  we 
obtain  the  following  source  encoding  intensity  trans¬ 
mittance  function. 


Mix, I)  v  it  !|X'' !  .  I  It') 

I  "  I 

when-  i \  ( I  >1  / 1„  is  tin-  diamelei  ol  a  i  iiciil.u  apcHmi1 

as  shown  in  Fig  3a, 

cjriMt  |i.  o-KI-.u- 

)  «•  i  |t).  otherwise. 
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/'  is  the  focal  length  of  the  collimating  lens  and  ).  is  the 
wavelength  of  the  extended  source.  As  a  numerical 
example,  we  assume  that  the  signal  si/e  is  /i„  -  5  mm, 
the  wavelength  is  A  =  5461  A,  focal  length  is 
/  =  300  mm,  then  the  diameter  iv  of  the  source  encod¬ 
ing  aperture  should  be  about  32.8  pm  or  smaller. 

We  now  consider  smeared  image  deblurring  [  i  3] 
operation  as  our  second  example.  We  note  that  the 
smeared  image  deblurring  is  a  I  -  I)  proces  :e.g  oper¬ 
ation  and  the  inverse  filtering  is  a  point-by-point 
processing  concept  such  that  the  operation  is  taking 
place  on  the  smearing  length  of  the  blurred  object. 

I  bus,  the  spatial  coherence  requirement  is  depending 
upon  the  smearing  length  of  the  blurred  object.  To 
obtain  a  physically  realizable  so  nice  encoding  func¬ 
tion,  we  let  the  spatial  coherence  function  al  the  input 
plane  P,  be 

n|.\, - \',|)  =  sine(  71  |.v,-.v',|).  (17) 

‘  -  'J.v2  *  'I 

where  ,1.v2  is  the  smearing  length.  A  sketch  of  (17)  is 
shown  in  Fig.  3b.  By  taking  the  Fourier  Transform  of 
(17),  we  obtain 

S(.vl)  =  rect|l'i,l|,  (18) 

where  tv  =  ( //)/(.  1.v2)  is  the  slit  width  of  the  source 
encoding  aperture,  as  shown  in  l  ig.  3b,  and 

rectf'-'U1’ 

|  »•  |  |(),  otherwise. 

For  a  numerical  illustration  if  the  smearing  length  is 
I  v_,  I  mm,  the  wavelength  is  >  5  lol  A.  and  the  focal 

length  is  /  300 nun.  then  the  slit  width  it  should  be 

about  161.8 pm  or  smaller. 

We  wouhl  now  consider  image  suhlraclion  |  14  |  lor 
our  third  illustration.  Since  the  image  subtiaetion  is  a 
I  I)  processing  operation  and  the  spatial  coherence 


requirement  i.-  depending  upon  the  corresponding 
point-pair  of  the  images,  thus  a  strictly  broad  spatial 
coherence  function  is  not  required.  In  other  words,  if 
one  can  maintain  the  spatial  coherence  between  the 
corresponding  image  points  to  he  subtracted,  then  the 
subtraction  operation  can  take  place  at  the  output 
image  plane.  Therefore,  instead  of  utilizing  a  strictly 
broad  coherence  function  over  the  input  plane  P,.  we 
would  use  a  point-pan  spatial  coherence  function. 
Again,  to  insure  a  physically  realizable  source¬ 
encoding  transmittance,  we  would  let  the  point-pair 
spatial  colicicncc  function  be  [10] 

/'<|v,  -a,|) 


where  2/i„  is  the  main  separation  of  the  two  input 
object  transparencies  at  plane  P2.  Np  I  a  positive 
integer,  and  we  note  that  »  <?(/.  liquation  ( P7|  repre¬ 
sents  a  sequence  of  narrow  pulses  which  occur  at 
jv2  —  v2l  =  n/i„.  where  n  is  a  positive  integer,  and  their 
peak  values  are  weighted  by  a  broader  sine  factor,  as 
shown  in  l  ie.  3c.  Thus,  we  see  that  a  high  degree  of 
spatial  coherence  is  ..laintamed  al  every  point-pair 
between  the  two  input  object  transparencies.  By  taking 
the  Fourier  transformation  of  tP.'T  we  obtain  the 
following  source  encoding  intensity  transmittance 

,V(|A|D-  V  reei||X|  "</| l ,  (20) 

vvhcic  u  is  the  slit  width,  and  <1  t/  / )  h„  is  the 
separation  between  the  slits.  It  is  clear  that  (201 
I'I'ii  .i  mI,  V  iimnliet  olTiaimw  slits  with  equal  spac¬ 
ing  (/,  as  shown  m  I  ig.  3c  As  a  numerical  example,  we¬ 
ld  the  sep, nation  of  the  input  objects  /t„  Idnim.  the 
wavelength  t  S46I  A.  the  focal  length  oi  the  col- 
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limaior  /  -  >00 min,  (liai  lilt  sp.il. nit  1/  between  the 
sills  is  Hi  4  pm  I  lit  slit  w  kill)  it  should  lit  smaller  lli.m 
il  2.  or  about  I  5  pill.  II  lilt  si/c  of  lilt  encoding  mask  is 
2  null  stpiait.  llitn  ilit  i  hi  in  lit  i  ol  sliis  \  is  ahoul  1 22 
I  lius  \\t  set  llial  w  illi  lilt  source  encoding  il  is  possible 
lo  intrtast  lilt  intensity  ol  die  illuiiiiuaiioii  \  lold.  and 
al  lilt  saint  tulle  it  maintains  lilt  point-pair  spatial 
coherence  requirement  lor  image  suhiratlion 
operation. 


Kxperiniental  Results 

In  this  section,  we  would  illustrate  (\u>  examples  as 
obtained  from  the  source  encoding  technique.  I  he  first 
experimental  illustration  is  the  result  obtained  for 
smeared  photographic  image  dcblurriug  with  encoded 
incoherent  source  as  shown  in  1-ig.  4  In  this  experi¬ 
ment  a  Xenon  art  lamp  vs  it li  a  green  interference  filler 
was  used  as  extended  incoherent  source.  A  single  slit 
mask  of  about  lOOptti  was  used  as  a  source  encoding 
mask.  I  lie  smeared  length  of  die  burred  image  was 
about  1  mm. 

figure  5  shows  an  experimental  result  obtained  from 
image  subtraction  operation  with  encoded  incoherent 
source.  In  this  experiment,  a  mercur>  art  lamp  with  a 
green  filter  was  used  as  an  extended  incoherent  source. 
A  mullislil  mask  was  used  to  encode  die  light  source. 
The  slit  width  iv  is  2.5  pm  and  die  spacing  between  sliis 
was  25  pm.  T  he  overall  si/.e  of  the  source  encoding 
mask  was  about  2.5  x  2.5  mm*.  The  mask  contains 
about  1(H)  sliis. 

from  these  experimental  results,  we  see  dial  the 
constraint  of  strictly  broad  spatial  coherence  require¬ 
ment  may  be  alleviated  with  source  encoding  tech¬ 
niques  so  that  il  allows  die  optical  information  pro¬ 
cessing  operation  can  he  carried  out  with  extended 
incoherent  source. 


mi 
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OPTICS 

b 

I  ii!  4.1  and  h  Pliologiaplnc  linage  dcblumng  with  encoded  exten¬ 
ded  incohcienl  xouice  i.n  Input  Mimed  object  and  tb)  deblurred 
linage 

('(inclusion 

We  have  derived  a  Fourier  transform  relationship 
between  the  spatial  coherence  function  and  the  source 
encoding  intensity  transmittance  function.  Since  the 
coherence  requirement  is  depending  upon  the  nature  of 
a  specific  information  processing  operation,  a  strict¬ 
ly  broad  coherence  requirement  may  not  be  needed 
in  practice.  The  basic  advantage  ol  the  source  encoding 
technique  is  to  alleviate  the  constraints  of  the  strict 
coherence  requirement  imposed  upon  the  optical  infor¬ 
mation  processing  system,  so  that  the  information 
processing  can  be  carried  out  with  encoded  extended 
incoherent  source.  I  lie  use  of  incoherent  source  tv' 
carry  out  the  optical  processing  operation  has  the 
advantage  of  suppressing  the  coherent  artifact  noise. 
In  addition,  the  incoherent  processing  system  is  usual¬ 
ly  simple  and  economical  to  operate.  Finally,  w  e  w  ould 
stress  that  the  source  encoding  technique  may  be 
extended  li>  white-light  optical  processing  operation,  a 
program  is  currently  under  investigation. 
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